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Abstract 



We propose a 3500 ton (3000 ton fiducial volume) Si02 neutrino detector with 
sampling calorimetry, charged particle tracking, and muon spectrometers to run in a 
Tevatron Fixed Target Program. Improvements to the Fermilab accelerator complex 
should allow substantial increases in the neutrino flux over the previous NuTeV quad 
triplet beamline. With 4 x 10^^ protons on target/year, a 5 year run would achieve 
event statistics more than 100 times higher than NuTeV. With 100 times the statistics 
of previous high energy neutrino experiments, the purely weak processes i^fi + e~ — > 

+ e~ and + e~ — > fg + fi~ (inverse muon decay) can be measured with high 
accuracy for the first time. The inverse muon decay process is independent of strong 
interaction effects and can be used to significantly improve the fiux normalization for 
all other processes. The high neutrino and antineutrino fiuxes also make new searches 
for lepton fiavor violation and neutral heavy leptons possible. In this document, we 
give a first look at the physics opportunities, detector and beam design, and calibration 
procedures. 
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1 Introduction 



The Neutrino Scattering on Glass (NuSOnG) experiment will consist of four detector 
modules, each composed of a finely segmented calorimeter followed by a muon spec- 
trometer. The detector will be illuminated by a neutrino or antineutrino beam from 
the Tevatron. In its five-year data acquisition period, NuSOnG will make precise mea- 
surements of three types of neutrino scattering and will accumulate the world's largest 
sample of electron-neutrino scatters. These data will provide unique opportunities to 
discover physics beyond the Standard Model (including, inter alia, lepton flavor vio- 
lation and new particles) as well as determine structure functions over a wide range 
of X and Q^. The breadth of anticipated measurements makes NuSOnG a program 
rather than an experiment; the design heritage ensures that the approach is low-risk 
and cost-effective. 

This Expression of Interest arises from our view that an experiment probing the 
high energy interactions of neutrinos is a necessary complement to the LHC and an 
important lead-in to the ILC. In the next few years, the LHC will reveal the nature 
of electroweak symmetry breaking; the Higgs mass will cease being a prediction of the 
electroweak theory and will become an input to the theory. Without the Higgs mass 
as a fit parameter, precision electroweak data, including neutrino scattering data, will 
be much more powerful as a tool for constraining that physics beyond the Standard 
Model which directly influences the electroweak sector. More important still, precision 
neutrino scattering will probe areas of phenomenology that may be inaccessible to the 
LHC and ILC. NuSOnG is not a precision test of the Standard Model; NuSOnG is a 
discovery experiment aimed at the terrain not covered by the collider experiments. 

This Expression of Interest presents the physics case and initial design for NuSOnG. 
The detector draws on the heritage of FMMF, CDHS, CHARM and CCFR/NuTeV. 
The design uses an Si02 target in one-quarter radiation length panels interleaved with 
active detector elements (proportional tubes and/or scintillator). This will provide the 
very high segmentation needed to ensure good separation between different classes of 
events. We will develop these ideas in the coming months and submit a proposal to 
the Fermilab Directorate. 

Our report is organized as follows: the physics opportunities follow in Section |2] 
Section [3] describes the flux and expected event rates; and Section |4] describes our 
preliminary design for the NuSOnG beam and apparatus. We summarize in Section |5j 
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2 Physics Opportunities 



The physics opportunities of the experiment arise from NuSOnG's uniquely high statis- 
tics: >20k neutrino-electron scatters and >100M neutrino-quark scatters. Roughly 
equal statistics will be obtained from antineutrino scattering. More information on the 
event rates for various processes is given in Sec. [3| These rates present a wide range 
of physics opportunities including precision electroweak measurements, direct searches 
for new physics, and parton distribution studies. 



2.1 Electroweak Precision Measurements 

NuSOnG's considerable discovery potential derives from its ability to do precision 
electroweak tests through two independent channels: electron scattering and quark 
scattering. These measurements probe for new particles and new neutrino properties 
beyond the present Standard Model. As examples, NuSOnG will be sensitive to extra 
Z bosons with masses beyond the 1 TeV scale (depending on the model), and to 
compositeness scales above 5 TeV. Thus the energy scales explored by this experiment 
overlap the LHC, and we present the discovery potential for the new physics we will 
explore within this context. This experiment also directly addresses questions raised 
by the "NuTeV anomaly," an electroweak precision measurement in disagreement with 
the Standard Model. 



2.1.1 Electroweak Measurements in Neutrino Scattering 



NuSOnG is sensitive to new physics through neutral current (NC) scattering. The 
exchange of the Z boson between the neutrino v and fermion / leads to the effective 
interaction: 



9l ^7m(1 - 75)1^ + Or i^7m(1 + 75)i^ 
5{/V(l-75)/ + 5^mi+75)/ 



where the Standard Model values of the couplings are: 

0, 

(4 - sin^ Ow 



9r 
9i 

9i 



or equivalently. 



9v 
9\ 



9l + 9r 



9l-9r 



(1) 



(2) 
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s'v = sI + s'r = v/p(;3'-2(?'sm'9w) , 



f f 

9l-9r 



V~P {4) 



(3) 



Here, /g and are the weak isospin and electromagnetic charge of fermion /, respec- 
tively. In these formulae, p is the relative coupling strength of the neutral to charged 
current interactions (p = 1 at tree level in the Standard Model). The weak mixing 
parameter, sin^ 6w, is related (at tree level) to to Gi?, Mz and a by 



sin^ 29w 



(4) 



NuSOnG is unique in its ability to test the NC couplings by studying scattering 
of neutrinos from both electrons and quarks. A deviation from the Standard Model 
predictions in both the electron and quark measurements will present a compelling 
case for new physics. 

Neutrino Electron Scattering 

The differential cross section for muon neutrino and antineutrino scattering from 
electrons, defined using the coupling constants described above, is: 



da 



2GlmeEn 



TT 



Wl9'v ± 9l9\f 



+{9l9'v T 9l9%f 



1 - 



dT 



-{{9l9vf - i9l9Af} 



rrifT dT 

~We: 



(5) 



The upper and lower signs corresponding to the neutrino and anti-neutrino cases, 
respectively. In this equation, E^, is the incident z/^ energy and T is the electron recoil 
kinetic energy. 

More often in the literature, the cross section is defined in terms of the parameters 
{9v^:9'a)^ which are defined as 



9!(r' ^ {Iglg-y) = p 
g^l ^ {2glg\) = p 

In terms of these parameters, we can write: 

dT 



2 sin 9w 



(6) 



da 



Gj^nieEu 
27r 



{9v ± 9Ar ^ + {9v T 9Tf 



dT 

En J En 



-{{97? -{97?} 



nieT dT 

~We: 



(7) 



When me <C Ei,, the third terms in these expressions can be neglected. If we introduce 
the variable y = T/Ey^ then 



da 
dy 



Gj^nieEn 
27r 



{9'(r' + 91^' + {9V - 97? {^ - y? 



(8) 
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Integrating over the region < y < 1, we obtain the total cross sections which are 



Note that 
Therefore, 



2-K 



1 



(9) 



{-I + Asm^ Ow) = /9^(l-2sin2^iy+4sin^0,y) , 



(10) 



a{u^ e) 
cj(P^ e) 



— P 

27r ^ 

27r 3 



16 

1—4 sin^ 9w + -IT sin^ 



1 -4sin2 6*^/ + 16sin^6'iy 



The ratio of the integrated cross sections for neutrino to antineutrino electron scat- 
tering is 



Re 



cr{i^fML e) 1-4 sin^ Ow + ^ sin^ 6w 

1 - 4sin2 9w + 16sin^ Ow 



Many systematics, including flux errors, cancel in this ratio, as does the p dependence. 
Fig. gtop) shows the results for sin 6w from many past experiments. 

NuSOnG will make independent measurements of the electroweak parameters for 
both and P^-electron scattering. We can achieve this via ratios or by direct extrac- 
tion of the cross section. In the case of f^-electron scattering, we will use the ratio of 
the number of events in neutrino-electron elastic scattering to inverse muon decay: 



NC 



X $^ 



P Ve 



a 



IMD 



X 



(13) 



Because the cross section for IMD events is well determined by the standard model, this 
ratio should have low errors and will isolate the EW parameters from NC scattering. In 
the case of running, the ratio is more complex because there is no equivalent process 
to inverse muon decay (since there are no positrons in the detector). In this case, we use 
the fact that, for low exchange energy in Deep Inelastic Scattering, the cross sections 
in neutrino and antineutrino scattering approach the same constant. A, as is explained 



in Sec. |3.3.2| Thus, for Deep Inelastic events with low energy transfer and hence low 

tIow Eh ad 



hadronic energy (5 < Ehad < 10 GeV), nI°^j^'^-^ = '^"A and A^"^^^""" = <^^A. The 
result is that we can extract the electroweak parameters to high precision using the 
ratio: 



KZs'''' , N{P,e- ^ i?,e-) _ <^ x '^^ 

X ' ~ — ~ ?\ 



P Ve 



a 



IMD 



X 



(14) 



The first ratio cancels the DIS cross section, leaving the energy-integrated v io v flux 
ratio. The IMD events in the denomenator of the second term cancel the integrated 
V flux. The NC elastic events cancel the integrated v flux. Alternatively, because we 



will have accurate knowledge of the flux as a function of the energy (see Sec. 3.3) we 
could directly measure the cross sections. 
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Figure 1: Measurements of siii^ 9w from past experiments. Top: neutrino-electron elastic 
scattering experiments. Bottom: neutrino DIS experiments. All DIS results are adjusted 
to the same charm mass (relevant for experiments not using P-W method). The Standard 
Model value, indicated by the line, is 0.2227. 
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An important point is that the two independent measurements, one in neutrino 
and the other in antineutrino mode, wih in turn ahow independent extraction of g'^ 
and gy'. The previous best measurement from and z^^ cross-section measurements 
is from CHARM II, which used 2677ib82 events in neutrino mode and 2752ib88 events 
in antineutrino mode [1] to find 

g'(p = -0.035 ± 0.012(stat) ±0.012(sys) (15) 
g-"/ = -0.503 ± 0.006(stat) ±0.016(sys). (16) 

This can be compared to electroweak measurements from LEP provide a very precise 
prediction of these parameters [2] : 

g!(,^ = -0.0397 ±0.0003 (17) 
g''/ = -0.5065 ±0.0001. (18) 

The CHARM II results are in agreement with LEP, but with large errors. Errors 
on the neutrino measurement must be substantially reduced in order to meaningfully 
probe for physics beyond the Standard Model. The goal of NuSOnG is to measure the 
neutrino-electron and antineutrino-electron cross sections to 0.7%. 



2.1.1.1 Neutrino Quark Scattering 

Substantially higher precision has been obtained using neutrino-quark scattering, 
which compares neutral-current (NC) to charged-current (CC) scattering to extract 
sin^ 9w- However, these experiments are subject to issues of modeling in the quark 
sector. Fig. [T](bottom) reviews the history of these measurements. 

The lowest systematic errors come from implementing a "Paschos-Wolfenstein style" 
[3] analysis, which would be the technique used by NuSOnG. This requires separated 
and i? beams, for which the following ratios could be formed: 

= ^ (19) 

R' = (20) 

(21) 

Paschos and Wolfenstein [3] recast these as: 

„iy „P nu „ pp 

where r = <y^ccl^^cc- ^ many systematics cancel to first order, including the 
effects of the quark and antiquark seas for u, d, s, and c. Charm production only enters 
through dvalence (which is Cabbibo suppressed) and at high x; thus the error from the 
charm mass is greatly reduced. The cross section ratios can be written in terms of the 
effective neutrino-quark coupling parameters gf^ and gj^ as 

= gl + rgl (23) 

R' = al + lal (24) 
R' = gl-gl = p\\-s\n^ew), (25) 
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Source 


Error 


Reduction in NuSOnG 


Statistics 


0.00135 


100 times the statistics 


z/e, z/g flux prediction 


0.00039 


see Sec. 


3.3.6| 


Interaction vertex position 


0.00030 


Better detector segmentation and 
more sophisticated shower identiflcation. 


1 j_l 11 

Siiower lengtii model 


0.00027 


Better segmentation and 
more sophisticated shower identiflcation. 


Counter efficiency and noise 


0.00023 


Better, Minos-style counter design 


Energy Measurement 


0.00018 


Likely to be at a similar level. 


Charm production, strange sea 


0.00047 


See Sec. 


2.4.4 




Rl 


0.00032 


Likely to be at a similar level. 




0.00022 


See Sec. 


2.4.5 




Higher Twist 


0.00014 


Likely to be at a similar level. 


Radiative Corrections 


0.00011 


Likely to be at a similar level. 


Charm Sea 


0.00010 


Under study 


Non-isoscalar target 


0.00005 


Glass is isoscalar 



Table 1: Source and value of NuTeV error on siri^Ow, and reason why the error will be 
reduced in the PW-style analysis of NuSOnG. 



in which 

gl ={2glglf + {2glgif= p\]^ - Ow + Isln"" Ow) (26) 

9r = {2glgl? + {2glgif = p\\sm'ew). (27) 

NuTeV fit for and B!^ simultaneously to extract sin^^iy, obtaining the value 
sin^ 9w = 0.2277 it 0.00162. The goal of NuSOnG is to improve on this error by a 
factor of two. Table [T] lists the errors which NuTeV identified and indicates those for 
which NuSOnG expects improvement. Many of the largest experimental systematics of 
NuTeV came from the method of separating CC and NC events, which relied on length. 
NuSOnG will have a more sophisticated model for differentiating CC and NC events, 
using shower shape and identification of Michel-electron followers from low energy pion 
decays. 

From Fig. [T] it is apparent that the NuTeV measurement is in agreement with 
past neutrino scattering results, although these have much larger errors. However, the 
NuTeV result is in disagreement with the global fits to the electroweak data which give 
a Standard Model value of sin^ 9w = 0.2227 [4j. Expressed in terms of the couplings, 
NuTeV measures: 

gl = 0.30005 ± 0.00137 (28) 
c/l = 0.03076 ±0.00110, (29) 

which can be compared to the Standard Model values of g\ = 0.3042 and g\ = 0.0301, 
respectively. Sec. |2.2| (below) considers possible sources for this disagreement, both 
within and outside the Standard Model. 
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2.1.2 NuSOnG and New Physics 



NuSOnG will provide important probes of physics beyond the Standard Model distinct 
from and complementary to those of the LHC. NuSOnG will seek indirect evidence for 
new physics by addressing anomalies in the precision electroweak data, and by pro- 
viding unique information about neutrino coupling to the Z. In addition, precision 
measurements from NuSOnG will help to disentangle the complicated set of observa- 
tions that will be present at the LHC and, in doing so, elucidate the mechanism of 
electroweak symmetry breaking. NuSOnG and the LHC provide distinct probes of 
new physics because new physics enters collider and neutrino scattering processes dif- 
ferently: neutrino physics measures different combinations of couplings to light quarks; 
neutrino scattering probes new physics at space-like momentum transfer (versus the 
time-like scattering at colliders); and systematics are very different between low and 
high energy experiments. Finally, NuSOnG will directly search for new particles and 
interactions in the lepton sector that might be missed by the LHC and must otherwise 
await discovery by the ILC. 

2.1.2.1 Nevj Physics Observed through Coupling to the Z 

NuSOnG is unique among experiments in its ability to address the nature of the 
neutrino couplings to the Z boson in the near future. In the Standard Model, the 
neutrino coupling to the Z- and PF-bosons is purely left-handed. Indeed, the fact 
that the neutrino coupling to the PF-boson and an electron is purely left-handed is, 
experimentally, a well-established fact (evidence includes precision measurements of 
pion and muon decay, nuclear processes, etc.). By contrast, the nature of the neutrino 
coupling to the Z boson is, experimentally, far from being precisely established [5]. 

The best measurement of the neutrino coupling to the Z-boson is provided by 
indirect measurements of the invisible Z-boson width at LEP. In units where the Stan- 
dard Model neutrino- Z-boson couplings are g'^ = 0.5, g'^ = 0, the LEP measurement 
[B] translates into (5^)^ + (5^)^ = 0.2487 ± 0.0010. Note that this result places no 
meaningful bound on g^^. 

Precise, model-independent information on g'^ can be obtained by combining I'^ + e 
scattering data from CHARM II and LEP and SLD data. Assuming model-independent 
couplings of the fermions to the Z-boson, + e scattering measures g'^ = 2p, while 
LEP and SLD measure the left and right-handed couplings of the electron to the Z. 
The CHARM II result translates into \gl\ = 0.502 ± 0.017 [5^, assuming that the 
charged-current weak interactions produce only left-handed neutrinos. In spite of the 
good precision of the CHARM II result (around 3.5%), a combination of all available 
data allows \gji/g1\ ~ 0.4 at the two a confidence level [5^. 

Significant improvement in our understanding of g'^ can only be obtained with more 
precise measurements of -|- e scattering, or with the advent of a new high intensity 
e~^e~ collider, such as the ILC. By combining ILC running at the Z-boson pole mass 
and at ^/s = 170 GeV, \g'^/gi\ ^ 0.3 could be constrained at the two a level after 
analyzing e~^e~ — > 7-l-missing energy events [5]. 

At NuSOnG, we estimate that g'^ can be measured at around the 0.86% level. This 
estimate is obtained by combining the statistical uncertainty (20,000 I' + e elastic scat- 
tering events) with an estimated 0.5% systematic uncertainty from the flux estimate. 
Fig. 12] (left) depicts an estimate of how precisely g'^ could be constrained if the Nu- 



ll 




Figure 2: Precision with which the right-handed neutrino-Z-boson couphng can be deter- 
mined by combining NuSOnG measurements of (/^ with the indirect determination of the 
invisible Z-boson width at LEP. In the left panel, we assume that the u + e scattering 
measurement is consistent with the Standard Model prediction g'^ = 0.5, while in the right 
panel we assume that the u + e scattering measurement is significantly lower, g'[ = 0.485, 
but still in agreement with the CHARM II measurement (at the one sigma level). Contours 
(black, red) are one and two sigma, respectively, while the star indicates the Standard Model 
expectation. See [5] for more details. 

SOnG result, assumed to agree with the Standard Model prediction, is combined with 
the indirect LEP constraints. One can clearly see that this measurement {\g^/g'i\ ^ 0.2 
at the two sigma level) compares favorably with the ILC capabilities described above. 
If the NuSOnG result is incompatible with Standard Model expectations but still in 
agreement with the CHARM II experiment, a combined NuSOnG-LEP analysis should 
be able to establish that ^ 0, as depicted in Fig. [2] (right). 

2.1.2.2 New Physics Observed through Oblique Corrections 

Precision neutrino scattering measurements made at NuSOnG can reveal new physics 
even when new particles are not created in the final state, through the effects of these 
particles in loops. For models of new physics in which the dominant loop corrections 
are vacuum polarization corrections to the gauge boson propagators ("oblique" cor- 
rections), the ST parameterization introduced by Peskin and Takeuchi [7] provides a 
convenient framework in which to describe the effects of the new physics. 

The ST parameterization begins with a reference Standard Model, including ref- 
erence values for the Higgs and top masses, and predictions for observables in this 
reference Standard Model. Differences between predicted and experimental values of 
the observables are then parameterized by and used to fit for S and T, which can then 
be compared to predictions from new physics. The full set of precision electroweak 
data can then be used to constrain 5 and T, as shown in Fig. |3j The T parameter is 
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S 

Figure 3: Three projected electroweak measurements from NuSOnG in S-T plane. LEP/SLD 
error ellipse is shown in red and the current NuTeV u — q measurement is shown as a light 
blue band. The ochre band shows NuSOnG z7 — e, the dark blue band shows NuSOnG 
u — q and the green shows NuSOnG u — e. The width of the bands correspond to 68% 
confidence level for statistics as described in the text. The NuSOnG measurements assume 
(5,T) = (0,0). 

sensitive to new physics that violates isospin and is zero for new physics that conserves 
isospin. Isospin-breaking new physics such as heavy non-degenerate fermion doublets 
or scalar multiplets would affect the T parameter. The S parameter is sensitive to 
isospin-conserving physics, such as heavy degenerate fermion doublets. 

The status of electroweak measurements are shown in Fig. [s] [H]. The combined 
analysis of the LEP and SLD data by the LEP Electroweak Working Group (EWWG) 
[5] indicates an allowed region shown by the small oval, centered at S" = 0.05 it 0.10 and 
T = 0.07 lb 0.11. A different choice of reference Higgs or top mass changes Standard 
Model predictions for observables and thus shifts the center of the ST plot [10] ; setting 
the Higgs mass to 1000 GeV would shift the center of the oval to roughly (5, T) = 
(0.12, —0.36). Measurements of the W mass, which are not shown, are also consistent 
with the LEP measurements. The highest precision neutrino result comes from vq and 
i>q scattering by the NuTeV experiment. This result clearly disagrees with the other 
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measurements, as discussed in Sec. 12.2 



The goal of NuSOnG is to make measurements which are competitive with or better 
than past electroweak measurements. These goals are indicated by the magenta ellipse 
and orange band on Fig. [3} The magenta ellipse shows the area in ST space which 
can be probed if a 0.7% measurement of the and z> NC electron-scattering cross 
sections is achieved. The orange band shows the improvement in the neutrino-quark, 
"Paschos-Wolfenstein" -style measurement which is expected from NuSOnG. 

Disregarding the NuTeV offset for the moment, one can now ask: how will this plot 
look in the era of LHC and what will NuSOnG add? We consider this question in light 
of three scenarios: 

1. a hght Higgs (115-200 GeV) 

2. a heavy Higgs (200-1000 GeV) 

3. no Higgs signal. 

2.1.2.3 NuSOnG Impact for a Light Higgs (115-200 GeV) Scenario 

A light Higgs is consistent with LEP/SLD and W mass data. The fit to the elec- 
troweak data excluding NuTeV indicates a mass less than 144 GeV at 95% CL. This is 
also consistent with the current best direct-search limit which finds mn > 114 GeV [9]. 
In the case of the lightest Higgs masses, where the cleanest signal may be in ^ 77, a 
clear observation above background will be experimentally difficult and may take some 
time. 

Once the LHC measurement of the Higgs mass is made, the center of the ST 
ellipse (Fig. [s]) will be fixed at a point (modulo any remaining uncertainty in the top 
mass). Our experiment is especially interesting if the NuSOnG result disagrees with 
this LEP-|-SLD-|-LHC point. If the LHC measurement is high, i.e. tuh ~ 200 GeV, the 
result would be marginally inconsistent with the M\y analysis, which is 85^28 GeV[9]. 
In this case, comparison with the scattering results from NuSOnG could resolve the 
question of a discrepancy between these measurements. 

If all other electroweak results are in good agreement, but disagree with NuSOnG, 
this would indicate new properties associated exclusively with the neutrino. An exam- 
ple would be decreased coupling of the neutrino to the Z boson, where suppression of 
the coupling comes from intergenerational mixing of the light neutrino with a moder- 
ately heavy neutrino: 

Uf, = (cos a)i^iight + (sin a)i^heavy (30) 
The Zv^v^ coupling is modified by cos^ a and the W^Vf^ coupling is modified by 



cos a. This model, inspired by the NuTeV anomaly (see Sec. 2.2), would yield a 
measurement in NuSOnG with a low NC-to-CC ratio in both the case of electron and 
quark scattering. 

These moderately heavy right-handed states, dubbed "neutrissimos" p[^, could 
have masses as low as just above the current bound of the Z mass. They may well be 
within the reach of the LHC and may appear as missing energy in events [12]. Some 
models allow for neutrissimos as light as ~ 100 GeV [13]. The neutrissimos decay very 
quickly, but not always invisibly. For example, in the reaction N ^ £ + W, the W 
may decay to either two jets or a neutrino-charged-lepton pair; only the latter case 
has missing energy. This may make recognition of the neutrissimo at LHC rather 
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difficult. In the case of mn < 130 GeV, a dominant decay mode of the Higgs (along 
with bb) could be into vN, where the neutrissimo subsequently decays. Reconstructing 
the Higgs in this case may be difficult at LHC; if neutrissimos exist, the result from 
NuSOnG may significantly improve our understanding of LHC results. 

With a large tuning among the neutrino Yukawa couplings [13j neutrissimos could 
be the seesaw right-handed neutrinos. Relatively "large" mixing is marginally consis- 
tent with other constraints, including neutrinoless double-beta decay, which constrains 
|[/e4p to be less than a few xlO~^ for a 100 GeV right-handed neutrino, and rare pion 
and tau decays, which constrain |C/^4|^ to be less than, most conservatively, 0.004 and 
|?7t4P to be less than 0.006. Other bounds come from ^ — > e conversion in nuclei and 
other charged-lepton-flavor violation. A new experiment to search for /i — > e has been 
proposed at Fermilab [H] should also be sensitive to neutrissimos. The combination 
of NuSOnG and this experiment will be powerful in identifying the existence of these 
particles. 

If the neutrissimo is a Majorana particle, it could be instrumental in elucidating 
the mechanism for leptogenesis. The present models of leptogenesis require very high 
mass scales for the neutral lepton, but theorists are pursuing ways to accommodate 
lower masses [15]. There also may be a wide mass spectrum for these particles, with 
one very heavy state required by standard leptogenesis models and others with masses 
in the range observable at LHC . 

2.1.2.4 NuSOnG Contribution in a Heavy Higgs (^200 — 1000 GeV) Scenario 

While present electroweak data excluding NuTeV favor a light Higgs (< 200 GeV), 
as indicated in Fig. [3] the Higgs mass can extend up to about 1000 GeV without 
violating unitarity [T7]. Thus, if LHC finds that the Higgs is between ~200 and 1000 
GeV and the LEP-I-SLD ellipse has no major systematic error, then new physics must 
explain the discrepancy. Candidate models of new physics may well affect the neutrino 
scattering and e~^e~ scattering differently, so the high-precision neutrino scattering 
measurements from NuSOnG will provide an important piece of the puzzle if the Higgs 
mass found at LHC is genuinely inconsistent with LEP-I-SLD predictions. 

Introduction of a fourth family would compensate for a modestly heavy (~ 300 
GeV) Higgs by shifting the LEP-I-SLD allowed region back up in S and T [18]. This 
family would need to exist above the bounds of direct searches, which is > 300 GeV. 
Mixing must be confined within the allowed bounds of the CKM matrix measurements 
[20\ . A nice feature of this model is that a fourth-generation Majorana neutrino could 
play the role of dark matter. Depending on the underlying physics, evidence of a 
fourth family would be apparent in a shift of the NuSOnG result on the ST plot. This 
could be especially important if the physics introducing the fourth family is from a 
mechanism like "Top See Saw" |21J, which will not be observable at LHC. The impact 
of this particular model on neutrino scattering is not yet thoroughly explored, but 
could prove interesting |22| . 

A classic method for masking a heavy Higgs is to introduce heavy Z bosons [23], 
which, as shown in ref. [lOj, tend to move the LEP-SLD ellipse upward in T, compen- 
sating for the heavy Higgs. Introduction of & Z tends to increase NC rate in neutrino 
scattering and also to move the neutrino result upward on the ST plot (although with 
a different dependence than the LEP-SLD result). 
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There are good theoretical reasons for considering the existence of additional neutral 
heavy gauge bosons. Extra Z bosons appear in various GUT and string-motivated 
extensions to the Standard Model [24] . For example, the -E(6) breakdown to 50(10) x 
results in the Z^. The 50(10) break down to S\J (5) x C/(l)x yields the Z^. Thus 
the new exchange boson could be: Z = Z^cos(3 + Z^sin(3, where the mixing angle 
P is an arbitrary parameter. Extra Z bosons also appear in other beyond Standard 
Model theories, including extra dimensions with gauge fields in the bulk [25]; little 
Higgs theories [25], which use heavy Zs to cancel divergences in the Higgs mass; and 
topcolor in which they drive electroweak symmetry breaking |27) . Heavy Zs provide 
a mechanism for new SUSY theories to evade the LEP bound of uih = 114 GeV |28] . 
These models all produce new physics signatures at LHC. The precision measurement 
from NuSOnG can aid in differentiating models. 

Models which introduce new physics to mask a heavy Higgs may seem contrived 
until one looks at the LEP+SLD data more closely. Up to this point we have considered 
the LEP+SLD measurements as a single result, however, many measurements enter 
this fit, and larger than expected inconsistencies between these measurements exist |29] . 
For example, there is a 3.2a discrepancy between the forward-backward {Afb) and left- 
right {Aiji) asymmetry measurements. Excluding the ApB result, the LEP+SLD fit 
yields mn < 115 GeV at 95%, with the best fit at 42 GeV - i.e. a range already 
excluded by direct searches, which require mn > 114 GeV at 95% CL. 

There are several ways to interpret this deviation. It may simply be that there 
are systematics involved in the ApB measurement which have yet to be identified and 
which would bring this result into agreement with the others. In this case, we are 
in the dramatic situation of having already ruled out the Higgs. The scenario of no 
Higgs is considered in the next section. Alternatively, new physics is involved. This 
result is dominated by purely leptonic measurements. On the other hand, the fit to 
the hadronic asymmetries, dominated by A'^p^ has two minima, at 450 and 3000 
GeV. Thus, one may either introduce new physics which produces a 20% shift on A^p^ 
alone; or introduce new physics which would indicate apparently low values of mn in 
the lepton-based measurements, when actually the value is large. Within any of these 
scenarios, new precision results from NuSOnG will be valuable for understanding the 
underlying physics. 



2.1.2.5 NuSOnG and the Case of No Higgs 



Higgsless models do not employ the Higgs mechanism to render the Standard Model 
renormalizable [30]; instead they introduce some other scheme. The Higgs mechanism 
enforces unitarity in the scattering amplitudes of longitudinally polarized gauge bosons, 
W^+Z1 W^Z^, for example. A requirement that the transition probability remains 
less than one gives the energy scale A at which a new mechanism must come into play, 

A ~ ^ l.STeV. (31) 

9 

Higgsless theories generally contain new mass bosons Vi with masses on the TeV scale 
that act to cancel the divergences in gauge boson scattering. Cancelling the ampli- 
tudes while respecting bounds from current electroweak couplings typically give small 
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couplings: 

VSM^Mw 



9WZV < = 0.04 (32) 



for Mf =700 GeV. 

At the LHC, the typical cross sections for Vi are hundreds of femtobarns, so, after 
cuts, the LHC experiments will record tens to hundreds of events in the first years 
of data taking. Since the Vi resonances serve the same purpose as the Higgs boson, 
additional information will be necessary to determine whether these resonances origi- 
nate from spontaneous symmetry breaking or from strong coupling between the known 
gauge bosons. The electroweak measurements from NuSOnG will play a role in under- 
standing the origin of such events, en route to a more complete explanation provided 
by the ILC. 



2.2 The NuTeV Anomaly 



The NuTeV anomaly is a 3a deviation of sin 9w from the Standard Model prediction 
[Ij. NuTeV employed the PW-inspired method discussed in Sec. 2.1.1.1 which resulted 
in a 0.75% measurement of the weak mixing angle (see Tab. [T|. Two systematic 
adjustments to the NuTeV result have been identified since the result was published. 
The first is the new measurement of the Ke3 branching ratio from KTeV, which does not 
significantly reduce the error, but introduces a correction moving the result away from 
the Standard Model. The second is the final measurement of the difference between 
the strange and antistrange seas (called "the strange sea asymmetry", see Sec. 2.4.4), 
which will pull the NuTeV result toward the Standard Model. A new analysis of the 
NuTeV data which will include these two corrections is expected be available in late 
summer, 2007 |32] . It should be noted that while an error from the strange sea appeared 
in the NuTeV analysis, no error on a strange sea asymmetry appeared in the original 
NuTeV analysis; this will be included in the upcoming re- analysis. 

NuTeV is one of a set of <C experiments measuring sin^ 9w. It was performed 
at Q2 = 1 to 140 GeV2, (Ql) = 26 GeV^, (Q?) = 15 GeV^, which is also the expected 
range for NuSOnG. Two other precision low measurements are from atomic parity 
violationf33] (APV), which samples ~ 0; and SLAG E158, a M0ller scattering 
experiment at average = 0.026 GeV^ Using the measurements at the Z-pole 
with = to fix the value of sin^ 9\y, and evolving to low Q^, Fig. |4] from ref. 
|31j . shows that APV and SLAG E158 are in agreement with the Standard Model. 
However, the radiative corrections to neutrino interactions allow sensitivity to high- 
mass particles which are complementary to the APV and M0ller-scattering corrections. 
Thus, these results may not be in conflict with NuTeV. The NuSOnG measurement 
will provide valuable additional information on this question. 

Since the NuTeV result was published, more than 300 papers have been written 
which cite this result. Various Beyond-the-Standard-Model explanations have been 
put forward; those which best explain the result require a follow-up experiment which 
probes the neutral weak couplings specifically with neutrinos, such as NuSOnG. Sev- 
eral "within-Standard-Model" explanations have also been put forward, based on the 
inherent issues involving scattering off quarks. NuSOnG can address these criticisms 
in two ways. First, we will provide better constraints of the quark-related distributions 
at issue. Second, we perform the measurement of the weak mixing angle in both a 
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purely leptonic mode (scattering from electrons) and via the PW method. Agreement 
between the two results would address the questions which have been raised. 



2.2.1 Explanations Within the Standard Model 

Four explanations for the NuTeV anomaly that are "within the Standard Model" have 
been proposed. These are: electromagnetic radiative corrections; higher order QCD 
corrections; isospin (or charge symmetry) violation; and the strange sea asymmetry. 
The radiative corrections will be disregarded here, since the results of this paper [36] 
are not reproducible. 

The effect of the possible explanations is illustrated in Fig. [5} On this plot, the 
solid horizontal line indicates the deviation of NuTeV from the Standard Model. The 
thick vertical lines, which emanate from the NuTeV deviation, show the range of pulls 
estimated for each explanation, as discussed below. The dashed horizontal line shows 
the estimated shift due the new i^es branching ratio. We do not yet have an estimated 
shift due to the new NuTeV strange sea measurement, but it is expected that this will 
move the dashed line toward the Standard Model [32]. 

Three "Standard Model" explanations may be considered next pj|l38]. First, the 
NuTeV analysis was not performed at a full NLO level; NuSOnG will need to undertake 
a full NLO analysis. But the effect of going to NLO on NuTeV can be estimated [39], 
and the expected pull is away from the Standard Model, as shown on Fig.|5j Second, the 
NuTeV analysis assumed isospin symmetry, that is, m(x)^ = d{x)^ and d{xY = u{x)^. 
Isospin violation can come about from a variety of sources and is interesting in its 



own right. NuSOnG's contribution to this study is discussed in Section 2.4.3 Various 
models for isospin violation have been studied and their pulls range from less than la 
away from the Standard Model to ~ la toward the Standard Model [3D]. We have 
chosen three examples [4D^ for illustration on Fig. |5j the full bag model, the meson 
cloud model, and the isospin QED model. These are mutually exclusive models, so 
only one of these can affect the NuTeV anomaly. Third, variations in the predicted 
strange sea asymmetry can either pull the result toward or away from the Standard 



Model expectation |4 H 1421 H3]. This issue is considered in detail in Sec. 2.4.4 



2.2.2 Beyond Standard Model Interpretations 

Chapter 14 of the APS Neutrino Study White Paper on Neutrino Theory [33] is dedi- 
cated to "the physics of NuTeV" and provides an excellent summary. The discussion 
presented here is drawn from this source. 

The NuTeV measurements of and , the NC-to-CC cross sections, are low. 
If one is assumes that the Higgs is light, then this must be interpreted as Beyond- 
Standard-Model physics that suppresses the NC rate with respect to the CC rate. Two 
types of models produce this effect and remain consistent with the other electroweak 
measurements: 1) models which affect only the Z couplings, e.g., the introduction of a 
heavy Z' boson which interferes with the Standard Model Z; or 2) models which affect 
only the neutrino couplings, e.g., the introduction of moderate mass neutral heavy 
leptons which mix with the neutrino. 



As discussed in Sec. 2.1.2.4, introduction of Z' bosons tend to increase the NC rate 
rather than suppress it. Thus there is only a small subset of models which produce the 
destructive interference needed to explain the NuTeV result. Models which introduce 
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Figure 5: Effect of various "Standard Model" explanations on the NuTeV anomaly. The 
y-axis is the deviation from the Standard Model. The solid hne is the NuTeV deviation. 
The dashed line is an estimate of the effect of correcting for the new Kg^ branching ratio. 
Thick black lines extending from the NuTcV deviation show the range of possible pulls from 
the various suggested sources, as described in the text. 
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a Z' which selectively suppresses neutrino scattering, without significantly affecting 
the other electroweak measurements, include cases where the Z' couples to i? — 
|45j or to — Lt [36]. In the former case, fitting the NuTeV anomaly requires that 
Mz' I gz' ~ 3 TeV. From the bounds from direct searches, this sets a limit on Mz' > 600 
GeV if the coupling is on the order of unity, but as low as 2 to 10 GeV if the coupling is 
~ 0.1%. The latter case is an example which improves the agreement between NuTeV 
and other results, but does not entirely address the problem. Its effectiveness in solving 
the NuTeV anomaly is limited by the data constraining lepton universality. This model 
addresses more than just the NuTeV anomaly. It is inspired by attempts to address 
bimaximal mixing in the neutrino sector. It has the nice features of also addressing 
the muon (g — 2) measurement and producing a distinctive dimuon signature at LHC. 

The case of models involving moderate-mass neutral heavy leptons, a.k.a. neutris- 
simos, have been discussed in the Sec. 2.1.2.3 and examples of viable models appear 
in ref. [11]. Eq. 30 described how the muon neutrino couplings might be modified by 
mixing. This idea can be extended to all three flavors, leading to a suppression factor 
for the Z coupling which is expressed as (1 — e^) and for the W hy {1 — e^/2), where 
£ = e, fj,, or r. This addresses the NuTeV anomaly and at the same time suppresses 
the invisible width of the Z, describing the LEP I data. 

If the NuTeV anomaly is due to Beyond Standard Model physics, then the effect 
will be visible in the neutrino-electron elastic scattering measurement also. Thus, if 
the NuTeV anomaly is borne out, NuSOnG would observe an ST plot similar to Fig.|6j 



2.3 Direct Searches for New Physics 
2.3.1 Light Neutrino Properties 

Evidence for three light neutrino masses has now been established through neutrino 
oscillations in solar, atmospheric, and reactor experiments (see references [U] through 
|61j). Furthermore, although the MiniBooNE experiment recently refuted the LSND 
two-neutrino oscillation scenario at Am? ~ 1 eV^ |62j, the question of the existence 
of multiple light sterile neutrinos still remains open [63] . These observations already 
require beyond-the-Standard-Model physics, and consequently raise phenomenological 
questions, such as: what are the mass and mixing parameters still allowed in sterile 
neutrino models? What do sterile neutrinos imply about neutrino mixing? Is the 
neutrino mixing matrix unitary, or is there effective freedom of mixing parameters? As 
we illustrate in the following sections, these are some of the questions that NuSOnG 
can potentially address. 



2.3.1.1 Matrix Freedom 



Perhaps the most interesting study of light neutrino properties which can be per- 
formed at NuSOnG is the search for evidence of "matrix freedom" or "nonunitarity." 
For example, in the case of existence of sterile neutrinos, the neutrino mixing matrix 
is extended to an x matrix, where A >3. Under that assumption, it has been 
suggested that the 3x3 part of the matrix describing the three active (SM) neutrinos is 
not necessarily unitary; or, equivalently, the three flavor eigenstates are non-orthogonal 
(the 3x3 neutrino mixing matrix is free) |64j- 
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Figure 6: Three projected electroweak measurements from NuSOnG in S-T plane for a 
model model with a heavy Higgs inspired by the NuTeV measurement [llj. In this model, 
{S,T) = (0.12, -0.36). The labeling is as in Fig. Isl 
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This introduces striking changes to the probabiUty formula for neutrino flavor tran- 
sitions. Assuming unitarity, the survival probability formula for a neutrino produced 
as flavor a is 

punAtary ^ ^ _ 4| [^^3 1 2 _ | 1 2] ^-^2 ^ (33) 

where one has made use of A31 = Am|;^^, and Am^i^ <C 1. In the case of matrix 
freedom, the mixing matrix is no longer unitary. The level at which unitarity is violated 
can be defined as Xa, where 



ajl^ = l-Xa, (34) 



with Xa being small. Under that assumption, the survival probability formula is then 
found to be 

pg^neral ^ punitary _ ^^^^^ _ ^^^^^2 ^-^2 ^^^j ^ ^2^ (35) 



As implied by Eq. 35 one of the main consequences of such scenario is instantaneous 
(-L =0) flavor transitions in a neutrino beam. This occurs regardless of the size of the 
mass splitting between the mostly sterile and mostly active states, and thus allows for 
a full-mass-range search for evidence of sterile neutrinos. A recent study [HS] suggests 
that current experimental data limit such an effect to up to the order of a few percent. 

As a result, several interesting and potentially observable phenomena can occur. 
Extending the argument of ref. [65], for instance, the non-orthogonality of and Ue 
that matrix freedom introduces, results in an instantaeous transition at L = from i/^ 
to fe ^5 . Thus one could observe an excess of events in a pure beam. 

The trick to searching for this instantaneous transition is to focus on an energy 
range where the Ve background is low and well constrained. In the case of NuSOnG, 
this is on the high energy tail of the flux, above E> 250 GeV. For the limits on 
transformation to fe[65j, which are at the ~ 1 x 10~^ level, NuSOnG would see an 
excess of ~ 200 z^g events in this high energy region. Fig. [7] shows the ratio of flux 
with transitions to flux without transitions. The abrupt cutoff is due to Monte 
Carlo statistics; higher energies can be explored. Assuming that such transitions indeed 
happen at the 10~^ level, one would expect up to a 10% increase in flux for E~ 350 
GeV. In that high energy region, the Ve flux is mainly from decay, which is well 
constrained by the i/^ events. Such an excess should therefore be measurable. 

Other interesting effects of matrix freedom [61] include the oscillatory behavior in 
the total (flavor-summed) CC event rate as a function of L/E, and (fake) CP- violating 
effects in the and u neutral-current event rates (the two rates oscillate differently 
with L/E). Potential observation of those effects at NuSOnG has not been explicitly 
considered at this stage, although it would be interesting to address this and we are 
planning to do so in the near future. Regardless of that, evidence of Ve contamination 
in a beam above expected background levels, something for which NuSOnG can 
search, would strongly support the matrix freedom hypothesis. 



2.3.1.2 Sterile Neutrino Oscillations 



Direct observation of sterile neutrino oscillations may also be possible in NuSOnG, 
depending on the mass and mixing parameters. Oscillations of active to light sterile 
neutrinos have been introduced to explain the LSND anomaly, as dark matter can- 
didates, and in describing the supernova collapse models. These ideas span a wide 
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Figure 7: Ratio of enhanced z/g flux due to transitions to u^. flux assuming no transitions. 
Obtained assuming lOOM z/^ deep inelastic scattering events. 



range of Am^ values. The LSND anomaly requires a sterile neutrino in the range 
of ~ 1 eV^ with moderate mixing (< 1%), while dark matter candidates and su- 
pernova collapse models require (> 1 keV)^. These models also require tiny mixing 
(10"^'^ < sin^ 26 < 10~^)[66 . NuSOnG probes an intermediate range of Am^, between 
the LSND and astrophysical allowed regions. However, since sterile neutrinos may 
come in families, it is worth exploring this previously uncharted territory. 

The NuSOnG experimental design consists of a 30-600 GeV muon neutrino beam, 
peaked at ~100 GeV, incident on a ~ 200-meter long detector located at L~1.5km from 
the neutrino source. This detector design allows for i^^ disappearance studies across the 
detector length by examining the scattering rate variation across the detector. Such 
searches would be limited by the detector energy resolution. Preliminary studies have 
shown that, assuming a 10% energy resolution. Am? ~ 600eV^ regions with mixing of 
< 0.1 can be probed easily. NuSOnG may also be able to explore smaller mixings and 
higher Am^s, depending on the final experimental design. 

NuSOnG can also probe for and z/g disappearance in the range of L/E = 
(1.5km/100 GeV) = 0.015, thus in the range of Am^ ~ 50 eV^. This is a range 
which has been covered by past experiments including CCFR |67,, CHDS [68 , and 
NOMAD |69]. However, the improved quality of the first principles prediction due 
to the new SPY secondary production data 
improvement of these limits. 



n^, discussed in sec. 3.3 should allow 



2.3.2 New Interactions 

2.3.2.1 Lepton Number Violation Searches 

The NuSOnG experiment possesses two valuable characteristics for the search for 
lepton number violation. First, it relies upon a high purity, high intensity beam as 
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its source of neutrinos; secondly, it employs an instrumented detector optimized to 
measure inverse muon decay with high accuracy. An experiment with these two features 
naturally lends itself to searches for the process: 



l^^l + e ^ n + z^e- (36) 

This interaction is forbidden by the Standard Model since it violates lepton family 
number conversation (ALg = —AL^ = 2). As such, observation of this reaction would 
immediately constitute direct observation of physics beyond the Standard Model. 

A number of theories beyond the Standard Model predict that lepton number is 
not a true conserved quantum number; this means that processes that violate lepton 
number are allowed to occur. Theories which incorporate multiplicative lepton number 
conservation [711 I72j . left- right symmetry f73], or the existence of bileptons [73] fall 
under this category. 

The differential cross-section for lepton-violating processes can be parametrized in 
the following form: 

^ = x9h^y . y(^y _r) + As-{l- r)), (37) 
ay IT 

where y is the fractional energy carried by the outgoing lepton, Gp the weak coupling 
constant, s the square of the center of mass energy of the system, and r the threshold 
factor, defined as rnj^/s. The parameters X,Av, and As describe the strength of the 
reaction and whether the process is vector or scalar in nature. It is typical to compare 
this process to that of inverse muon decay: 



\.{Av{—-^)+As). (38) 



The signature for such a reaction is the tagging of an during antineutrino run- 
ning with the same signature as expected from inverse muon decays. The main back- 
grounds to this reaction include (a) contamination, (b) contamination, and (c) 
charge misidentification of candidate events. Our current estimates place a very small 
beam contamination during antineutrino running: about 0.4% contamination of v^s 



and a 2.3% contamination of Ue and T>(, neutrinos (See Sec. 3.1). Charge misidentifica- 
tion is expected to be very small, on the order of 10~^. If we assume a conservative 
knowledge of the backgrounds at the 5% level, this would imply a limit on the lepton 
number violation cross-section ratio of better than 0.2% (at 90% C.L.) for V-A cou- 
plings and less than 0.06% for scalar couplings. Previous searches, based on 1.6 x 10^^ 
protons on target and smaller target masses, have placed limits on this cross-section 
ratio to less than 1.7% at 90% C.L. for V-A couplings and less than 0.6% for scalar 
couplings [75] . The NuSOnG experiment can therefore reach an improvement of over 
an order of magnitude compared to previous searches. This limit can be improved 
if further selection criteria are used in removing unwanted beam impurities or the 
quasi-elastic background contamination. 

2.3.2.2 Inverse Muon Decay 

The study of inverse muon decay, z/^ -|- e~ n~ + fg provides access to the helicity 
structure of the weak interaction distinct from muon decay experiments. The weak 
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interaction polarizes the incident f^, making inverse muon decay an excellent place to 
study departures from V — A couplings. For inverse muon decay, a cx {g^''^g^''^)'^{l — e) 
|76| where e = /i — (— 1) and h is the helicity of the incident muon. Ref. [77] has 
measured e < 4.1 x 10~^ and the current limit on = {g^^'^g^'^) > 0.96 [TH]. For 
a measurement of the total cross section scaled to the predicted cross section, the 
uncertainty on the coupling is g^i = {l/2)aa-/crsM- 

For NuSOnG, we expect > 200k inverse muon decay events, which would give a 
statistical uncertainty of 0.002 on g^i^- However, we will need to determine the neutrino 
flux. Taking the + e~ ^ + cross section as known gives the neutrino flux to 
0.7%. Since we plan to use the inverse muon decay events for determining the flux for 
the electroweak measurements, NuSOnG will need to measure the efficiency and fiducial 
volume for both processes to better than 0.7%. Combined with other systematics, we 
should be able to achieve an total uncertainty of about 1-2% on gYi^ ^.n improvement 
by a factor of four. 

The key background will come from CCQE events that have small hadronic energy. 
We expect our high granularity will allow us to keep the systematic error from this 
source well below 1%, but this needs study. 

Obviously, the manner of analysis described above is somewhat questionable. Ul- 
timately, one would want to carry out a combined analysis of both neutrino elastic 
scattering on electrons and quarks and of inverse muon decay in the context of a spe- 
cific model which relates the charged and neutral current coupling constants. For such 
an analysis, 1-2% uncertainty should still be achievable. 

2.3.3 New Particles 

2.3.3.1 Long-lived, Light Neutral Heavy Leptons 

Another interesting NuTeV result arose from the search for long-lived, light (< 15 
GeV) neutral heavy leptons. This was performed in a helium-filled decay region located 
upstream of the calorimeter. In the mass region of 2.2-15 GeV, NuTeV has a small 
expected background (0.07 ± 0.01 events), but observed three events. All events had 
two muons originating from a vertex within the helium decay region and missing energy. 

Since publication in 2001, no widely accepted explanation has been found. In 2006, 
DO published a search for a similar decay signature in proton-antiproton interactions 
|80j . No events were found and some production models were excluded. The most 
viable remaining model is by Dedes et al, which hypothesizes that the events are from 
decay of long-lived neutralinos. These are produced in the NuTeV beam dump through 
B hadron decays [81,. No other experiment has been able to match NuTeV's running 
conditions to further explore this intriguing result. 

NuSOnG can address the question by including a low-mass (helium-filled) decay 
region between the calorimeter segments. Assuming parameters similar to those of 
NuTeV (except for a 20-fold increase in the number of protons on target), NuSOnG 
would expect to see 60 events with an expected background of 1-2 events. The sen- 
sitivity would scale directly with the decay volume, so the increased length compared 
to NuTeV (26 m ~40 m) would increase this to 90 signal events over a 2-3 event 
background. Observing no signal would finally settle this outstanding question. 
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These decay regions allow exploration for a signal from a beyond-the-Standard- 
Model particle in other decay modes as well; other interesting modes include fin, fie, 
en and ee. NuSOnG's sensitivity to other new particles is similarly improved over 
NuTeV by the increase in beam intensity and decay volume, allowing us to study new 
regions of phase space. 

2.3.3.2 Muomc Photons 

In the mid-1990's there was interest in searching for "leptonic photons" - mass- 
less vector particles that couple according to flavor. Electronic, muonic, and tauonic 
photons, 7e, 7;^, and 7,- were introduced f82]. Production occurs in secondary meson 
decays such as tt ^ 1/^/27^, and detection can proceed through 7^ + Z ^ /i"*" + + Z, 
where Z is the charged nucleus. These events have small missing px compared to the 
"trident" background, v + Z ^ v + ^"^z^" + Z. The search by CHARM II sets the best 
limit at 1.6 X 10"^ [83] . 

Since this time, neutrino oscillations have been confirmed (see references [37] through 
[61j). This complicates the theory of "muonic photons," since, in this case, lepton 
flavor-charge is not conserved. As pointed out in reference |82j, a theory with a non- 
conserved charge cannot have massless vector particles and a Coulomb-like potential. 
It appears very difficult to evade this problem. 

Nevertheless, NuSOnG should search for these events. With higher rate and better 
segmentation than CHARM II, NuSOnG should have sensitivity in the range of ~ 10^''. 
A significant excess would be quite startling. 

2.4 IVIeasurement of Parton Distribution Functions 

The Deeply Inelastic Scattering (DIS) process provides crucial information about the 
structure of the proton which is used to determine the Parton Distribution Functions 
(PDFs). For example, in the recent CTEQ6HQ analysis, DIS data accounted for more 
than two-thirds of the data points used in the analysis]^ As such, the DIS measurements 
form the foundation for the many calculations which make use of the PDFs. 

In the basic DIS process, leptons scatter from hadrons via the exchange of an 
intermediate vector boson: {7, , Z}. Different boson probes couple to the hadrons 
with different factors, and it is important to combine data from these different probes 
to separate the different flavor components in the hadron. Unfortunately, three of the 
four DIS probes {W^,Z} have a (relatively) large mass and couple only weakly; this 
introduces a number of complications: 

• The statistics for these weak processes are limited as compared with the photon- 
exchange processes. 

• To compensate for the weak cross section, typically heavy nuclear targets {e.g., Fe 
and Pb) are used; this introduces nuclear corrections when the results are scaled 
from the heavy target back to proton or isoscalar targets. 

The NuSOnG experiment will generate high statistics (> lOOM DIS events) mea- 
surements on an intermediate atomic- weight nuclear target (Si02). This will provide 

^Specifically, there were 1333 DIS data points used out of the 1925 total. [84] 
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precise information on the linear combinations of PDFs which couple to the weak 
charged currents (W^), which can significantly improve the parton distribution fits. 
In this section, we first introduce the basics of DIS and the connection to parton distri- 
bution functions. Then we concentrate on three aspects of parton distribution studies 
where NuSOnG can make a unique contribution to the physics: 

• Improved understanding of nuclear effects in neutrino scattering. 

• Study of Charge Symmetry Violation 

• Measurement of the Strange Sea 

• Measurement of a'^ and a'^ 

The latter two items are directly relevant to the electroweak studies proposed for 



NuSOnG (see Sec. 2.2.1). 



2.4.1 Deep Inelastic Scattering and Parton Distribution Functions 

The differential cross section for neutrino DIS depends on three structure functions: 
F2, xF^ and R^. It is given by: 



dxdy 
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(39) 



where the it is +(— ) for yiv) scattering. In this equation, x is the Bjorken scaling 
variable, y the inelasticity, and the squared four-momentum transfer. 

The function xF^{x, Q^) is unique to the DIS cross section for the weak interaction. 
It originates from the parity-violating term in the product of the leptonic and hadronic 
tensors. For an isoscalar target, in the quark-parton model. 
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Defining xF^ 



l{xF^^ + xF^^), at leading order in QCD, 



xFs^LO = ^ xq{x, 

i=u,d.. 



xq{x, 



(40) 
(41) 

(42) 



To the level that the sea quark distributions have the same x dependence, and thus 
cancel, xFs can be thought of as probing the valence quark distributions. The difference 
between the neutrino and antineutrino parity violating structure functions, A(a;F3) = 



xF^ — xF^ , probes the strange and charm seas. 

Analogous functions for F2{x,Q'^) and Rl{x,Q'^) appear in both the cross section 
for charged lepton (e or fi) DIS and the cross section for v DIS. At leading order. 



F2,LO= Yl e\xq{x,Q^)+xq{x,Q^)), 



(43) 



i=u,d.. 
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where e is the charge associated with the interaction. In the weak interaction, this 
charge is unity. For charged-lepton scattering mediated by a virtual photon, the frac- 
tional electromagnetic charge of each quark flavor enters. Thus F!^^ and F^^^^^ are 
analogous but not identical and comparison yields useful information about specific 
parton distributions [87J. Rl{x, Q^) is the longitudinal to transverse virtual boson ab- 
sorption cross-section ratio. The best measurements for this come from charged lepton 
scattering rather than neutrino scattering. In the past, neutrino experiments have used 
the charged lepton fits to Rl as an input to the measurements of xF^ and ^2^85]. This, 
however, is just a matter of the statistics needed for a global fit to all of the unknown 
structure functions in x and bins |86]. With the high statistics of NuSOnG, precise 
measurement of Rl will be possible from neutrino scattering for the first time. 

In addition to fitting to the inclusive DIS sample, neutrino scattering can also probe 
parton distributions through exclusive samples. A unique and important case is the 
measurement of the strange sea through opposite sign dimuon production. When the 
neutrino interacts with an s or d quark, it produces a charm quark that fragments into 
a charmed hadron. The charmed hadron's semileptonic decay (with branching ratio 
Be ~ 10% ) produces a second muon of opposite sign from the first: 

ly^, + N — > fi- + c + X (44) 
^ s + n+ + u^. (45) 

Similarly, with antineutrinos, the interaction is with an s or d, 

z7^ + N^//++ c + X (46) 

fi- + Uf,. (47) 

The opposite sign of the two muons can be determined for those events where both 
muons reach the toroid spectrometer. Study of these events as a function of the kine- 
matic variables allows extraction of the strange sea, the charm quark mass, the charmed 
particle branching ratio {Be), and the Cabibbo-Kobayashi-Maskaka matrix element, 

\Vcd\. 

For a more in-depth review of precision measurement of parton distributions in 
neutrino scattering, see ref. [88] . 

2.4.2 Nuclear Effects 

Historically, neutrino experiments have played a major role in expanding our under- 
standing of parton distribution functions through high statistics experiments such as 
CCFR [85J, NuTeV [89], and CHORUS [90^. However, the high statistics extract a 
price since the large event samples require the use of nuclear targets - iron in the case 
of both CCFR and NuTeV and lead in the case of the Chorus experiment. The problem 
is that if one wants to extract information on nucleon PDFs, then the effects of the 
nuclear targets must first be removed. NuSOnG can provide key measurements which 
will improve these corrections. 

In the case of charged lepton deep inelastic scattering, there are data available from 
nuclear targets covering the range from deuterium through iron and beyond. Thus, 
it has been possible to perform detailed studies of the j4-dependence as a function of 
X and from both the cross section and the structure function F2. Such is not the case 
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in V and V interactions where the corrections can be different for both cross sections 
or, equivalently, for F2 and xF-^. In this case one must rely on theoretical models 
of the nuclear corrections. This is an unsatisfactory situation since one is essentially 
measuring quantities sensitive to the convolution of the the desired PDFs and unknown 
- or model dependent - nuclear corrections. 

It is important to address the question of nuclear effects in neutrino scattering so 
that the neutrino data can be used in fits without bringing in substantial uncertainties. 
For example, in a recent analysis [HI] the impact of new neutrino data on global fits for 
PDFs was assessed. The conclusion reached in this analysis was that the uncertainties 
associated with nuclear corrections precluded using the neutrino data to constrain the 
nucleon PDFs. If this uncertainty is addressed, the neutrino data will be a powerful 
addition to these fits. 

Furthermore, nuclear effects are interesting in their own right. Comparison of 
the charged and neutral lepton scattering data can provide clues to the sources of 
the major features which appear in nuclear effects: shadowing, antishadowing, and 
the EMC effect. There is phenomenological evidence which suggests that the nuclear 
corrections for the v and v cross sections might be rather similar and, in both cases, 
somewhat smaller than the corresponding corrections in charged lepton deep inelastic 
scattering. These latter two observations differ from the pattern suggested by the 
theoretical model |92] for nuclear corrections used in the analysis. 

Fig. [sjshows some results from Ref. [911 in the form of "data/theory" averaged over 

and presented versus x. The results are from a global fit but are plotted without 
the model-dependent nuclear corrections which were used in the fits. What is striking 
is the similarity of the v and V results, and the overall pattern of deviations, similar 
to that seen in charged lepton DIS, although the deviations from unity are somewhat 
smaller. It is interesting to note that there is no clear indication of the turnover at low 
X which is observed in charged lepton scattering, called shadowing. However, this may 
be due to kinematic limits of the measurements. 

To make progress in understanding nuclear corrections in neutrino interactions, 
access to high-statistics data on a variety of nuclear targets will be essential. This will 
allow the j4-dependence to be studied as a function of both x and Q^, as has been done 
in charged lepton deep inelastic scattering. PDFs from global fits without the neutrino 
data can then be used to make predictions to be compared with the ^-dependent 
1/ and z7 cross sections, thereby allowing the nuclear corrections to be mapped out for 
comparison with theoretical models. 

The primary target of NuSOnG will be Si02. However, we can address this issue 
by replacing a few slabs of glass with alternative target materials: C, Al, Fe, and 
Pb. This range of nuclear targets would both extend the results of Miner z^a to the 
NuSOnG kinematic region, and provide a check (via the Fe target) against the NuTeV 
measurement. 

Given the NuSOnG neutrino flux, we anticipate 58A; i^-induced and 30/c z^-induced 
CC DIS events per ton of material. A single ton would be sufficient to extract F2{x) 
and xF^^x) averaged over all Q^; a single 5 mx5 mx2.54 cm slab of any of the above 
materials will weigh more than that. The use of additional slabs would permit further 
extraction of the structure functions into separate (x, Q^) bins as was done in the 
NuTeV analysis, at the potential expense of complicating the shower energy resolution 
in the sub-detectors containing the alternative targets; this issue will be studied via 
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Figure 8: Comparison between the reference fit and the unshifted Chorus and NuTeV 
trino data without any nuclear corrections. 
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Material 


Mass of 
2.54 cm slab (tons) 


Number of slabs needed 
for NuTeV-equivalent statistics 


C 


1.6 


33 


Al 


1.9 


27 


Fe 


5.5 


10 


Pb 


7.9 


7 



Table 2: Alternative target materials for cross-section analysis 



simulation. 

Table |2] shows that two 50-module stacks would be sufficient to accumulate enough 
statistics on alternative nuclear targets for a full structure-function extraction for each 
material. However, for basic cross-section ratios in x, a single slab of each would suffice. 

2.4.3 Isospin Violations 

When we relate DIS measurements from heavy targets such as 26-^^ 82'^-'^^ back 
to a proton or isoscalar target, we generally make use of isospin symmetry where we 
assume that the proton and neutron PDFs can be related via a u <^ d interchange. 
While isospin symmetry is elegant and well-motivated, the validity of this exact charge 
symmetry must ultimately be established by experimental measurement. There have 
been a number of studies investigating isospin symmetry violation \93\ IMl [95l [96] : 
therefore, it is important to be aware of the magnitude of potential violations of isospin 
symmetry and the consequences on the extracted PDF components. For example, the 
naive parton model relations are modified if we have a violation of exact p <-> n isospin- 
symmetry, (or charge symmetry); e.g., Un{x) ^ dp{x) and Up{x) ^ dn{x). 

Combinations of structure functions can be particularly sensitive to isospin viola- 
tions, and NuSOnG is well suited to measure some of these observables. For example, 
residual u, ^-contributions to AxF^ = xF^^ — xFJ^ from charge symmetry violation 
(CSV) would be amplified due to enhanced valence components {uy(x),dv{x)}, and 
because the d ^ u transitions are not subject to slow-rescaling corrections which 
strongly suppress the s ^ c contribution to AxFs. [95j Here the ability of NuSOnG to 
separately measure xF!^ and xF^ over a broad kinematic range will provide powerful 
constraints on the sensitive structure function combination AxFs. 

There are a wide variety of models that study CSV [931 EH ES |96]. One method 
to quantify possible CSV contributions is via a one-parameter "toy" model where we 
characterize the CSV as a rotation in isospin space: q^^^ = Nq Y^^, Rqqi{6) q'p, where R 
is a rotation matrix, and Nq is the normalization factor. For example, the u-distribution 
in the neutron can be related to the proton distributions via: 

u^s^(x, g2) = [cos2(0) upix, Q2) + sin2(6') dp{x, Q^)] . (48) 

For 9 = 7r/2, we recover the symmetric limit Up{x, Q^) = dn{x, Q^). While this param- 
eterization does not offer any explanation for the source of the CSV, it does provide 
a simple one-parameter (9) characterization which is flexible enough to quantify the 
range of CSV effects. (For more details, cf. Ref. [95].) 

At present, there are constraints on isospin violation from a number of experiments 
which cover different ranges of x and Q^. For example, we note that while the above 
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"toy" model leaves the neutron singlet combination {q + q) invariant at the < 2% level 
in the region x e [0.01; 0.1], it would lower the NC observable [^{u + u) + ^{d + d)^ in 
this region by about 10%. An effect of this size would definitely be visible in the NMC 
F2/F2 data which has an uncertainty of order a few percent .jSTj The measurement 
of the lepton charge asymmetry in W decays from the Tevatron |97| |98] places tight 
constraints on the up and down quark distributions in the range 0.007 < x < 0.24. 
While only strictly telling us about parton distributions in the proton, these data 
rule out isospin violations at the 5% level, as demonstrated in Ref. In addition, 
there are also fixed-target Drell-Yan experiments such as NA51 [99] and E866 |100] 
which precisely measure d/u in the range 0.04 < x < 0.27; these are also sensitive to 
isospin- violating effects. 

NuSOnG will be able to provide high statistics DIS measurements across a wide x 
range. Because the target material (Si02) is very nearly isoscalar, this will essentially 
allow a direct extraction of the isoscalar observables. Consequently, if isospin violations 
are present, they can be measured more precisely than would be the case on a highly 
non- isoscalar target. 



2.4.4 Measurement of the Strange Sea 

There are several reasons why an improved measurement of the strange sea is of inter- 
est. First, it contributes to the low-Q^ properties of the nucleon in the nonperturbative 
regime - a question of practical as well as intellectual interest, since many precision 
oscillation experiments are being performed in the 1 to 20 GeV (hence, nonperturba- 
tive) range. It is critical for charm production which provides an important testing 
ground for NLO QCD calculations. In addition, understanding the threshold behavior 
associated with the heavy charm mass is of interest to future neutrino experiments. 
Distinguishing the difference between the s{x) and s{x) distributions, 

xs~{x) = xs{x) — xs{x), (49) 

is even more important, and poses additional challenges. First, it is of intrinsic interest 
in nucleon structure models ^43l |38l |371 IIOH 1102] . Second, the integrated strange sea 
asymmetry, 

= s-{x)dx, (50) 







has important implications for the precision measurement of the weak mixing angle 



in deep inelastic scattering of neutrinos (c/. Sec. 2.2 and references IH] US] [ST] |38l 
SH lion US] ) • This was not recognized at the time of the NuTeV sin^ 9w publication; 
an error due to S~ will be included in the NuTeV reanalysis, to be presented in late 
summer 2007 [52]. 

Historically, information on the s(x) and s{x) distributions was derived from in- 
clusive cross sections for neutral and charged current DIS via A^xF^). These analyses 
made the implicit assumption that the s{x) and s{x) seas had the same distribution 
in X. Because the strange sea is relatively small compared to the dominant u{x) and 
d{x) processes, the resulting uncertainties on the strange sea components were large. 
For example, the strangeness content of the nucleon, as measured by the momentum 
fraction carried by s or s, is of order 3% at Q = 1.5 GeV. For this reason, the strange 
PDF was typically parametrized using the ansatz s{x) = s{x) = k{u + d)/2, where n 
measured the deviation from SU{2>) fiavor symmetry at some low value of Q. 
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Figure 9: xs^{x) vs x aX Q"^ = 16 GeV^. Outer band is combined errors, inner band is 
without Be uncertainty. 

Introducing information from opposite sign dimuon production allows substantial 
improvement in the strange PDF measurement. Neutrino-induced dimuon produc- 
tion, {v/i>)N fj,'^ij,~X, proceeds primarily through the sub-processes W~^s — > c and 
W~s — > c (respectively), so this provides a mechanism to directly probe the s(x) and 
s{x) distributions without being overwhelmed by the larger u{x) and d{x) distribu- 
tions. Hence, the recent high-statistics dimuon measurements |103| I104| I105| I106| 1107] 
play an essential role in constraining the strange component of the proton. 

The highest precision study of s~ to date is from the NuTeV experiment |108] . 
The sign selected beam allowed measurement of the strange and antistrange seas inde- 
pendently, recording 5163 neutrino-induced dimuons, and 1380 antineutrino-induced 
dimuon events in its iron target. Figure [9] shows the measured asymmetry between 
the strange and antistrange seas. With more than 100 times the statistics of NuTeV, 
NuSOnG will have substantially finer binning. 

The integrated strange sea asymmetry from NuTeV has a positive central value: 
0.00196±0. 00046 (stat) ±0.00045 (syst) Iq^ooioi (external). The "external" error on the 
measurement is dominated by the error on the average charm semi-muonic branching 
ratio. Be which is determined by other experiments. This error currently is about 10%. 
A rescan of Chorus data, which would increase the statistics, is under consideration 

m- 

The key to an improved result on the strange sea from NuSOnG is in a more precise 
measurement of Be at NuSOnG energies. This can be accomplished in two ways. First, 
the very high statistics of NuSOnG allow for an accurate fit to B^ and the s and s 
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Figure 10: World measurements of B^. 

distributions simultaneously. Second, we plan to incorporate an emulsion detector 
into the design. The emulsion will be scanned by the Nagoya University group. This 
group has substantial expertise, having provided the emulsion and scanning for Chorus, 
DoNuT and other emulsion-based experiments. The goal will be to obtain > 10k events 
in the emulsion during the NuSOnG run. 

Beyond this, we will also consider placing a liquid argon TPC of similar size to mi- 
croBooNE |llUj (70 tons fiducial volume) or even Gargamelle (20 tons fiducial volume) 
in the gap between two of the NuSOnG subdetectors to directly measure Be- If one 
were to, for example, use a volume comparable to the Gargamelle bubble chamber, we 
could observe on the order of one million charged current events within it for 5 x 10^^ 
POT. This would yield approximately 100,000 events with charm in the final state, 
and about 10,000 dimuon events. 

In addition to an improved measurement of Be, the more finely-grained liquid argon 
TPC and/or emulsion detectors could be used to aid the calibration of the four glass 
detector modules by measuring any differences between hadron and electron showers 
from pion and electron beams versus those within a neutrino induced event. Coupled 
with the precision test beam, it may also be possible to improve understanding of the 
background due to muons produced by pion and kaon decays in the hadron shower. An 
improved parameterization of this background, currently from a CCFR measurement 
|11H I112| I113j could help extend the kinematic range of charmed dimuon measurements 
beyond what was possible for the NuTeV and CCFR experiments. 

2.4.5 Measurement of the Total Cross Section 

Precision measurement of the total neutrino and antineutrino cross sections at high 
energies will be valuable to a future neutrino factory experiment which seeks to make 
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Figure 11: World measurements of the total u and u cross sections. See references [Hn] and 
[TTi] through [T^ . 



precision measurement of CP violation. Because NuSOnG can measure the flux to 
0.5% (see Sec. 3.3, precise measurements can be made. Also, higher accuracy on the 
ratio of ja^ will also improve the electroweak measurement (see Tabjl]). 

Fig. 11 shows (t/ Ey for the muon neutrino and antineutrino charged-current total 
cross-section as a function of neutrino energy ( |89| and |114j - P^H] l. The error bars 
include both statistical and systematic errors. The results are from a wide range 
of target materials, but the experiments with the smallest errors and largest energy 
range used iron. The straight lines are the isoscalar-corrected total cross-section values 
averaged over 30-200 GeV as measured by the experiments in Refs. |115j to |117j . The 
flt [129] gives: a^'^'^/Ey = (0.677 ± 0.014) x l{)-^^cw? /GeV] a^^'°/E^ = (0.334 ± 
0.008) X 10~^^ cm^/ GeV. The average ratio of the antineutrino to neutrino cross- 
section in the energy range 30-200 GeV is ^-^ ^^"^^i/ /so _ 0.504 it 0.003 as measured by 
Refs. [89] and |114j - [TT7] . Note the change in the energy scale at 30 GeV. 

The most precise measurements are systematics limited. The largest contributions 
to the systematics in recent experiments (CCFR, NuTeV) come from flux normaliza- 
tion, the model parameterization used in determination of the flux, and the charm mass 
used to parameterize charm threshold. NuSOnG measures the neutrino flux normaliza- 
tion to high precision via the IMD events. Also, as discussed in Sec. 3.3.2 NuSOnG's 
high statistics allow cuts which substantially improve the model parameterization er- 
ror. Lastly, the charm mass, rUc is expected to be improved from the high statistics flts 
to the opposite sign dimuon events described in the previous section. While more study 
is needed, it likely that NuSOnG can substantially improve on the world measurements 
of the total cross section and the cross section ratios. 
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Figure 12: NuSOnG flux in neutrino mode (left) and antineutrino mode (right). Black: 
muon neutrino flux, red: muon antineutrino flux, blue: electron neutrino and antineutrino 
flux 

3 Neutrino Flux and Event Rates 



3.1 The Neutrino Flux 

For the purposes of this expression of interest, we assume the same SSQT design as 



was used at NuTeV. The resulting neutrino (antineutrino) flux [33] is shown in Fig. 12 



left (right). The flux is shown in black, z^^ in red, and Ue + i^e in blue. The shape of 
the flux is dominated by the dichromatic neutrino spectrum from vr and K two-body 
decay. 

In neutrino mode, 98.2% of neutrino interactions are due to tt"'" and secondaries, 
while in antineutrino mode 97.3% come from vr" and K~ . The "wrong sign" content 
is very low, with an 0.03% antineutrino contamination in neutrino mode and 0.4% 
neutrino contamination in antineutrino mode. The electron-flavor content is 1.8% in 
neutrino mode and 2.3% in antineutrino mode. The major source of these neutrinos is 
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iTgg decay, representing 1.7% of the total flux in neutrino mode, and 1.6% in antineu- 
trino mode. Other contributions come from K^eS, KseS, charmed meson, muon, Ac, 
A, and S decays. 

Precise knowledge of the electron-flavor content is crucial for many NuSOnG anal- 
yses. The largest source of error in the knowledge of the electron-flavor content in 
NuTeV was from the branching ratio, which led to an error on Ue content of 1.4% 
|33j . While the other sources of z^gS have large fractional errors, they constitute a much 
smaller fraction of the flux. An error of 1.5% for the electron-flavor contamination, 
consistent with NuTeV, will be assumed for NuSOnG. 



3.2 Event Rates 

The approximate event rates presented here serve to set the scale for the physics case 
presented in this document. They are based on running the Nuance event generator 
1 130] with the NuTeV flux, and then scaling to the expectations of NuSOnG assuming 
a 3 kton fiducial mass. Some simplifying assumptions, which will be corrected as the 
simulation develops, have been made. For example, C2 is used as a target rather than 
Si02. Also, note that Nuance is not yet tuned as a high energy event generator. Thus, 
these event rates are only representative. 

For neutrino running, approximate event rates for 5 x 10^^ protons are: 
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21k neutrino — electron NC elastic scatters 

OA; neutrino — electron CC quasielastic scatters {IMD) 

The above were run for and beams. The relative ratios of event- weighted 
contents in neutrino mode are: - 98.33%, - 0.08%, - 1.56%, Ve 0.03%. The 
relative ratios of event-weighted contents in antineutrino mode are: - 0.42%, i/^ - 
98.07%, i/e - 0.26%, Ve 1.26%. 

3.3 Precision Measurement of the Flux from Events in 
the Detector 

Precise knowledge of the neutrino flux is key to many of the physics goals of the 
experiment. The goal, which is ambitious, will be to measure the neutrino flux as a 
function of energy to a precision better than 0.5%. This goal is a design-driver for the 
experiment. In this section, we outline an analysis plan to achieve this goal using the 
event types described in the previous section. 

The flux will be determined through the following steps: 

1. The inverse muon decay (IMD) events (f^ -|- e~ — > /i~ -|- Ve) are, in principle, 
ideal for measuring the total flux because the IMD cross section is well known in 
the Standard Model. Therefore, these events will be used to determine the nor- 
malization of the flux. An important background to this measurement, however, 
comes from the CCQE events (z^^ -|- n — > ii + p), which must be subtracted. In 
this step, the predicted number of CCQE events is based on external cross section 
measurements. The error on the external cross section is likely to be the limiting 
systematic on the normalization determined in this step. 

2. The shape of the flux is measured using the traditional "fixed measurement 
method, which was applied in CCFR [85. 1114] . and NuTeV [89^, and is currently 
being used for in the Minos Experiment |131j to measure the shape of both the 
neutrino and antineutrino flux. The flux shape is then normalized by the IMD 
events from step 1 to obtain the initial flux prediction. 

3. The initial flux prediction is used to determine a more precise CCQE cross section 
based on the NuSOnG data. 

4. Step 1 is repeated using the more precise cross section determined in step 3. 
This produces the final normalization which is used to scale the results of step 2, 
yielding the final flux. 

When the analysis is performed, it may be more effective and efficient to combine the 
above steps into a single multiparameter fit to the IMD and CCQE data, constrained 
by the external cross section information. However, for transparency we will consider 
the stepwise approach below. 

Reaching the goal of < 0.5% systematic error depends mainly on the systematics 
of the IMD total event rate measurement. To set the scale of the problem, the best 
measurement of IMD events to date, from CHARM II, had a systematic error of 3% 
|132] . Thus, we must achieve an order of magnitude improvement in the IMD total 
systematic error. While we present a well-grounded back-of-the-envelope argument 
below, this level of measurement has yet to be demonstrated by simulation. That is a 
priority for future work on the development of NuSOnG. 
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Figure 13: Momentum distribution of protons in u CCQE events, from the Nuance Event 
generator. 

Two useful cross checks of the measured NuSOnG flux are possible. First, one 
can extract an energy binned neutrino flux from the IMD events in step 1. Because 
of angular resolution, this flux may have substantial smearing, but it can be used as 
a compelling cross-check of the flux shape derived in step 2. Second, the neutrino 
to antineutrino flux ratio can be compared to the first principles prediction based on 
secondary production measurements. 

3.3.1 Step 1: The IMD Measurement for Normalization 

NuSOnG expects to observe > 200k IMD events during neutrino running. The high 
statistics is a consequence of both the high neutrino flux and high neutrino energy. High 

energy is required because the threshold for IMD scattering is Ei, > > 2^ = 10.9 
GeV. The SSQT beam design for NuSOnG produces minimal flux below 30 GeV, well 
within the range of IMD production. This indicates that there will be high statistics 
for IMD events in all flux bins. 

These events will be used for total flux normalization, with the shape determined 
using the Fixed v method described in step 2. This is done because, while these events 
can in principle be fully reconstructed assuming that the incoming neutrino enters 
parallel to the z-axis, the reconstruction in practice suffers substantial smearing. At 
~100 GeV, IMD events will have scattering angles of < 1 mrad. This is similar in 
magnitude to the expected divergence of the beam, which was 0.62 mrad in NuTeV. 
Angular resolution errors are expected to be at a similar level. 

IMD events must be separated from background, mainly due to CCQE-like interac- 
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tions (which include both real CCQE interactions and single vr events where the pion 
was absorbed in the nucleus, and thus are effectively CCQE events). IMD events are 
qualitatively different from CCQE-like ones in two ways: there is no hadronic energy 
in the event, and there is a strict kinematic limit on the transverse momentum of the 
outgoing muon, px < 2mgE^. It is therefore crucial to design NuSOnG for observation 
of very low hadronic energy in the presence of a muon track, and for excellent angu- 
lar resolution on the outgoing muon. The fine segmentation of NuSOnG should allow 
hadron identification in the presence of a muon to substantially lower energy (~0.75 
GeV) compared to 1.5 GeV for Charm II |132] and 3 GeV for NuTeV[75j. One can 
see from the momentum distribution of the protons produced in CCQE interactions, 
shown in Fig. 13 that this approach will allow NuSOnG to cut more of the CCQE 
background than was possible in the previous experiments. NuSOnG also expects bet- 
ter IMD resolution than Charm II, due to the finer segmentation, which will reduce 
backgrounds. 

Events which produce very low energy pions can also produce a background, al- 
though at a lower level than the CCQE background. NuSOnG's open trigger will allow 
many of the low multiplicity DIS events and CCvr^ events to be identified and cut due 
to the presence of subsequent michel electrons which come from the tt~^ /x"*" — > e"*" 
decay chain. A 50 MeV michel electron will traverse 12 cm of glass, producing hits in 
up to four chambers in the vicinity of the interaction vertex. 

The IMD method for determining the flux proceeds in the following manner. After 
cutting on hadronic energy, minimum energy for the outgoing muon, and no michel 
electrons near the vertex, the plot of muon pT will show a sharp peak at ~ 
superimposed on a broad continuum of background events extending to high px- The 
continuum is fit and extrapolated under the IMD peak, to extract the number of IMD 
events. This is divided by the theoretical cross section to yield the flux. 

At the high energies of NuSOnG, the only nuclear effect expected for CCQE events 
comes from the Pauli exclusion effect. This produces an overall suppression of the cross 
section across all energies. Both the NuTeV and Charm II measurements suffered from 
the lack of availability of precise information on the Pauli exclusion effect. This resulted 
in an error on the Charm II measurement from the CCQE model of 2.1%. 

NuSOnG will be in the fortunate position that a number of new measurements of the 
CCQE cross section on nuclear targets will be available as inputs into the CCQE model. 
Results from MiniBooNE |133j and SciBooNE |134j will address Pauli suppression in 
CCQE interactions on carbon. Minerva [135] is studying a series of nuclear targets, 
and are willing to consider running a glass target for NuSOnG, if we were to supply 
the target panels. The precision on the CCQE cross section in the NuSOnG era may 
be 5%, which is ~5 times better than the CHARM II era measurements. Thus, at this 
step, the CCQE model error for NuSOnG may be as low as 0.4%. 

A CCQE model error which was not addressed in the Charm II analysis was the 
long-standing discrepancy between models and data at low |133] . Low events 
having small scattering angles represents a significant error on the extrapolation under 
the IMD peak. This discrepancy has recently been resolved by MiniBooNE under a 
dipole form- factor model [133] . Minerva plans to address the dependence of the 
form factor in a model-independent way |136] . We will assume that the discrepancy 
will be fully addressed by the time of the NuSOnG run. 

Another 1.5% systematic error in Charm II came from the model of the other sources 
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of low hadronic energy CC events, which are dominated by A resonant production. 
As described above, NuSOnG expects a substantially lower contamination from these 
sources because of the lower energy threshold and the michel electron veto. For those 
background events which are not cut, the modeling of these sources will substantially 
improved using Minerva data. In Minerva, the tracks from CCvr"'" events are well 
reconstructed, so models pT distribution of outgoing muons can be tuned. Similarly, 
Minerva offers the opportunity to accurately parametrize CCpi'' events. Again, the 
total error on the low hadron multiplicity events in 2015 is expected to be on the order 
of 5%, so the modeling of these backgrounds should not be a limiting systematic. 

IMD events can be cut from the sample due to electromagnetic radiation of the 
muon near the vertex region which is mistakenly identified as hadronic energy at the 
vertex. The NuTeV IMD analysis [75] assigned a 1% systematic error due to radiative 
effects. We address this error in two ways. First, in the NuTeV experiment, the photons 
immediately converted in the 10 cm iron plates, while in NuSOnG, the 2.5 cm glass 
plates are 0.25Ao, giving photons a 50% probability of traversing three plates before 
showering. Fewer IMD radiative events will therefore be misidentified as events with 
hadronic energy at the vertex. Second, because of the higher segmentation, NuSOnG 
will employ an improved model of electromagnetic showers, reducing the systematic 
error. 



3.3.2 Step 2: The Fixed-z/ Measurement to Determine the Shape 

The central premise of the Fixed v method for measuring the flux is that, for small 
hadronic energy exchange (u), the differential cross section is independent of energy 
to a good approximation. The Fixed-z^ method utilizes this fact to measure the rela- 
tive flux between energy bins and the relative flux between neutrino and antineutrino 
interactions. External input is then needed to determine the overall normalization. 

To motivate the premise, consider the differential cross section at a fixed v inte- 
grated over all x: 

da ..^ B V C v"^ . 

^(l + T^-TTT^)- (51) 



In this equation. 



A = k J F2{x,Q^)dx, (52) 
B = -k j[F2{x,Q^)TxF3{x,Q^)]dx, (53) 



l+2Mx 



Mx 



C = B-k F2{x,Q'){ - l)dx, (54) 

J l + R{x,Q'') V 

(55) 

where k = {G'^pM)/-K^ and =p refers to neutrinos (— ) or antineutrinos (-|-). For sim- 
plicity, first consider i/ — > 0. The cross section becomes equivalent to A, which is a 
constant. Since it is impossible to measure scattering for = 0, consider scattering 
for V = vq where vq <^ E^. As long as is small enough, the terms which depend on 
v^jEy will have negligible contribution. Thus for a fixed, low value of da/du A, 
independent of beam energy. Note that terms B and C differ for neutrinos and an- 
tineutrinos. However, as long as vq/Ei, is negligible, these terms do not contribute and 
the cross section for antineutrinos is equal to the cross section for neutrinos. 
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From this, one can see how to measure the relative fluxes. If one measures the 
number of events at a given uq in bins of E^, one can solve for the flux: 

^E,) = N{E,,iJo)/A . (56) 

The relative change of flux between two energy bins is independent of A: 

<!>{E'^^^)/<!>{Elf''^) = N^'''\E^,uo)/N{E'^'^'^,uo) . (57) 

Since the neutrino and antineutrino cross sections are equal, this method also allows 
the relative fluxes to be extracted, independent of A. 

^E,)/^Ep) = N{E,,uo)/N{E^,uo) . (58) 

Thus one can extract the relative bin-to-bin and neutrino-to-antineutrino fluxes strictly 
from the data, with no theoretical input on the value of A. 

In practice one uses a low region, deflned hy v < vq where vq is some appropriate 
upper limit. CCFR and NuTeV, used v < uo = 20 GeV, which allowed high statistical 
precision for the measurement. From the theoretical point of view, however, this was 
not optimal since the goal was to measure the flux down to E^, = 30 GeV, thus at 
V = 20 GeV, the v/Ey terms were not negligible. The flux is then given by: 

/.^'0 dN_ 

^iE.)= —^±^^du . (59) 

A flt to dN/du determines B /A and C /A. One can test the quality of the bin-to-bin 
result by fltting a /E to a. line. A good fit results in small slope, due to QCD effects on 
the order of a few percent (somewhat smaller in antineutrino mode), with small error. 
NuTeV found values consistent with expectation [89]: 



A(f) 
A(f) 



(-2.2±0.8)%/100GeV, (60) 
(-0.2±0.8)%/100GeV. (61) 



The NuTeV analysis indicated a good fit to a straight line, as expected. The extracted 
shape of the flux was obtained to very high precision across the full energy range by 
this approach. 

NuSOnG has an important advantage over NuTeV when implementing this method, 
in that the high statistics and good segmentation will all us to reduce this range of 
the low v substantially, perhaps to as low as u < vq = 10 GeV. This should allow an 
even more precise measure of the shape than was obtained by past experiments, since 
the contribution of the flt to the B and C terms will be reduced. In particular, the 
systematic error contribution from the charm mass will be substantially reduced. 

NuTeV also required v > b GeV to cut the resonance region. NuSOnG is also likely 
to introduce such a cut. However, this should be revisited in light of the expected new 
data from Minerva in the resonance region. 

The most important detector systematic to this measurement is likely to be the 
muon energy scale. NuTeV achieved knowledge of the muon energy scale to 0.7%, 
although the absolute calibration beam was known to 0.3%. The difficulty was mapping 
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across the full area of the toroids. For NuSOnG to achieve its goal of measuring the 
flux with < 0.5% total error, the muon energy scale will need to be known to about 
0.25%. Careful thought must be put in to understand how to achieve this. 

In past experiments, the next step was to obtain the absolute flux by normalizing 
to the world's total, which is a/E^ = 0.667 ± 0.014 x 10-38cmVGeV. It necessarily 
introduces a 2% normalization error into this method. NuSOnG will use the IMD 
events to perform the absolute normalization, rather than relying on the world average 
neutrino cross section measurement. This is done by scaling the total flux measured 
in neutrino mode with the Fixed v method to Y^j^N^^^{Ei) J a^^^dE. At the end 
of this step, the predicted flux is expected to be known to ~ 1%. 

3.3.3 Step 3: A Precise Measurement of the CCQE Cross Section 

At this point in the procedure, the limiting systematic is likely to be the CCQE cross 
section model error in the IMD normalization. In this step, this cross section is further 
constrained using the CCQE data in NuSOnG. 

The background to the CCQE cross section analysis will be the low hadronic energy 
events. These can be reduced using the michel veto method discussed in Step 1. Beyond 
this, because CCQE scatters extend to higher angles, excess hits due to the presence 
of charged pions and photons from tt*^ decay should be more easily resolved from the 
photon track. NuSOnG expects ~ 500k CCQE events, and thus stringent cuts can 
be applied to remove backgrounds without substantial statistical error, assuming the 
efficiency of the cuts can be well-understood. 

The goal will be for NuSOnG to measure the CCQE cross section to <2%. This 
would be a very valuable measurement in its own right, as well as allowing for im- 
provement in the flux extraction in the following steps. This result can be used to 
constrain the normalization for a glass-target measurement in Minerva. Ratios to the 
other nuclear target cross section measurements by Minerva then allow precisely deter- 
mined measurements at low E across a wide range of nuclei. This will be useful input 
to future precision neutrino oscillation measurements. 

3.3.4 Step 4: The Final NuSOnG Flux 

Once the CCQE cross section has been determined at the <2% level, one can iterate the 
IMD analysis of step 2 and then renormalize the distributions in step 3. The resulting 
flux is expected to have errors of < 0.5%. 

3.3.5 Cross Checks 

Two useful cross checks of the flux are possible. The first takes the measured flux and 
compares it to the IMD event rate in energy bins. The second uses external data to 
cross check the shape and normalization of the antineutrino flux. 

The first cross check compares the shape of the neutrino flux determined at step 1 
to the shape determined through step 2. This will be done by running the final flux 
through the MC and using it to predict the IMD rate in energy bins. We will then 
extract the predicted flux in energy bins to compared to the measurement performed 
in step 1. This provides a powerful consistency check. 
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We can also cross check the fluxes obtained by the above method using a first- 
principles prediction based on external secondary production measurements. The ab- 
solute predictions in neutrino and antineutrino mode are unlikely to be an effective 
cross check because of large errors in the secondary production predictions, which vary 
from 5 to 10%. However, the prediction of the ratio of the neutrino to antineutrino 
fluxes may be possible to high precision. This requires some investigation. 

Reference |139j provides a compendium of secondary production experiments in 
Table 3. None extend up to 800 GeV. The most relevant experiment was NA56/SPY 
at 450 GeV, which took data on beryllium targets [70]. This experiment published 
yields of tt^ ,tt~ ,K~^ and K~ with errors on each measurement of ~ 5%. However, 
because many of the systematics cancel in ratio, the tt~ /pi~^, K~ /K^ and tt/K ratios 
are each determined to ~ 2.5%. This data should allow a good cross check of the 
individual vr and K shape contributions. We may choose to run for a short period at 
450 GeV in order to have an exact cross-comparison. 

3.3.6 The Electron Neutrino Flux 

We will begin by tuning the NuSOnG Beam Monte Carlo using the recent secondary 
meson production measurements described above. The new K production results will 
improve the first principles prediction for electron neutrinos beyond those of NuTeV. 
The electron neutrino contamination then can be further constrained through the pre- 
cision measurement of the flux, which can be tied to the Ue flux, and through the 
measurement of v^. CCQE events. 

Once the muon neutrino flux is measured to high precision, it can be used to 
constrain the electron neutrino flux. This is because the Ve (i^e) background is largely 
due to {K~) decays in neutrino (antineutrino mode). Using the measured peak 
from events, the beam Monte Carlo can be precisely tuned. Having measured the 
CCQE cross section precisely in the process of determining the u^^ flux, this result can 
then be applied to Vf, CCQE events to cross check the Ve flux prediction. 
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4 Preliminary Design 



This report focuses upon the determination of the physics goals of the experiment. In 
order to maintain reahstic goals, we have developed a preliminary design for a beam and 
detector based on existing technology. There are two particularly challenging aspects 



of the design. The first is the high Tevatron intensity discussed in sec. 4.1 The second 



is the high precision required for the detector calibration discussed in sec. 4.3.3 



The 2007 Fermilab Steering Group Report considers the Tevatron-based neutrino 
beam described here. The preliminary concept for the facility recieved an endorsement 

m- 



4.1 Proton Delivery to NuSOnG 

Our goal is to obtain 2 x lO^'^ protons on target during a 5-year run. This section 
outlines how we might achieve this goal. 

Proton delivery occurs via the following lines: 

• The Linac 

• The Booster 

• The Main Injector 

• The Tevatron 

• Extraction to targeting 

The existing Linac and the Booster should perform to the level needed by NuSOnG 
without problems. The Booster fills the MI in batches of 5 x 10^^ protons and will 
operate between 9 and 15 Hz by 2015. The Proton Plan projects 7 x 10^^ protons in 
each MI fill by 2010 |138j . Two pulses from the MI are used to fill the Tevatron. In 
principle, therefore, it is conceivable that the Tevatron could receive nearly 1.5 x 10^^ 
protons per fill under this scenario. 

Let's suppose that with care the Tevatron can accelerate 8 x 10^^ ppp to 800 GeV 
using two pulses from the Main Injector at 4 x 10^'^ each pulse, similar to today's MI 
operation. To date, the highest intensities extracted from the Tevatron in a single pulse 
at 800 GeV were around 2.5 to 3 x 10^^. The limiting issue was longitudinal instabilities 
for energies above 600 GeV at high intensities, as the bunch length shrank. "Bunch 
spreaders" were used to compensate. A better method to compensate will be required 
for NuSOnG. However, advances in rf techniques and technology and in damper systems 
make finding a satisfactory solution conceivable. More detailed study is needed. 

Our proposal is for a Tevatron cycle time of 40 s, with a 1 s flattop at 800 GeV. 
Since the MI cycle time will be 2.2 s, and we need two injections, our impact on NuMI 
is 4.4/40 = 11% of their run time. 

If the uptime for the Tevatron is 66%, then we will receive 5 x 10^ cycles per year. 
At 8 X 10^^ PPP, this gives 4 x 10^^ protons per year. We then achieve our goal in five 
years of running. 
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4.2 Neutrino Beam Design 



4.2.1 Target 

Beryllium oxide was the target material in NuTeV and prior high energy neutrino 
beamlines |140| I141j . Beryllium is efficient at producing secondary mesons, and BeO 
has good structural and thermal properties. The NuTeV target consisted of two 30 cm 
long, 2.5 cm diameter segmented rods in a copper cooling block, mounted on a movable 
drive that could select between centering the beam on either of the two targets or no 
target. This target was designed to accept up to 1x10^'^ protons per pulse (ppp). A 
similar target will be acceptable for NuSOnG, but it may be a challenge to provide 
adequate cooling at our design intensity of 8 x 10^^ protons per cycle. The NuTeV 
protons were delivered in five 4 msec "pings" separated by 0.5 sec; we intend to have 
one pulse of about 200 msec. This means that our instantaneous heating rate will be 
somewhat lower than NuTeV 's, but the total number of protons per cycle is eight times 
higher. In NuTeV the beam width was 0.6 mm, which was significantly smaller than 
necessary; a wider more diffuse beam would help relieve the localized heating problem. 
Careful design of the target support and cooling system will be a necessity. 

4.2.2 SSQT 

A Sign Selecting Quadrupole Train (SSQT) can be used to provide beams of either 
neutrinos or antineutrinos with very low contamination from either wrong-sign muon 
neutrinos or electron neutrinos from neutral kaons. The NuTeV SSQT utilized two 
dipoles and six quadrupoles, with two dumps |140| I141j . The first dipole provided 
a 6.1 mrad bend for 250 GeV daughter mesons of the selected sign. In antineutrino 
mode the unreacted protons are bent in the opposite direction and are absorbed in 
the first dump. In neutrino mode the protons are absorbed in the second dump. The 
first two quadrupoles capture the secondary beam. A second dipole then bends the 
beam by another 1.6 mrad, enhancing the sign separation and sweeping out low energy 
particles produced by scraping in upstream magnets. Neutral particles are not bent and 
therefore travel away from the detector. NuSOnG will use a similar SSQT. The only 
challenge will be designing the proton dumps for our significantly higher intensity. In 
the NuTeV upstream dump in antineutrino mode, the dump temperature approached 
100 C at 1.3 X 10^^ protons per pulse; the temperature limit was 110 C. The NuSOnG 
dumps will need to be water-cooled. 

4.2.3 Monitoring 

Primary beam monitoring in NuTeV was accomplished with four beam position mon- 
itors (BPMs), four vacuum segmented wire ionization chambers (SWICs), four sec- 
ondary emission electron detectors (SEEDs), a beam current toroid, and a thin foil 
secondary emission monitor (SEM) [TIOl [MI] . The toroid, SEM, SEEDs and BPMs 
measured proton intensity; the BPMs, SWICs, and SEEDs monitored position. It was 
found that the SEM degraded over the course of the run, so the beam toroid was used 
as the primary measure of intensity. The BPMs and SEEDs gave closely correllated 
position measurements, and the SWICs and SEEDs gave beam profiles that agreed well 
except in the tails; the SEED tails dropped more rapidly than those from the SWICs. 
With the exception of the SEMs, which would suffer even more radiation damage at 
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Parameter 



Value 



Total target mass 3.492 

Fiducial mass 2.975 kt 

Total length 192 m 

Number of glass planes 2500 

Number of toroid washers 96 

Number of muon detector wire planes 60 
(two coordinates each) 



Table 3: Summary of NuSOnG detector parameters. 



our higher intensities, a combination of any of these monitoring devices could be used 
by NuSOnG. 

4.3 Detector Design 

This section details our first ideas about the detector configuration; these are summa- 
rized in Tab. |3j In thinking about NuSOnG, we have drawn on previous large, high 
energy neutrino detectors whose characteristics are summarized in Table |4} NuSOnG 
represents a natural evolution of these designs and we believe this makes construction 
low risk. Of particular note regarding Table |4] is the excellent performance achieved 
by CHARM II using digital proportional tubes, a glass target, and fine granularity. 
The primary event signatures NuSOnG will need to identify are: 

• charged current deep inelastic scattering, characterized by a hadronic shower and 
a high energy muon 

• neutral current deep inelastic scattering, characterized by a hadronic shower 

• inverse muon decay, + e~ — > /u~ +i'e, which is characterized by a high energy 
muon accompanied by no hadronic activity. 

• neutrino and antineutrino electron scattering, characterized by an electromagnetic 
shower with no hadronic activity 

• stopped muon decay, which results in an electromagnetic shower with energy up 
to 50 MeV. These events will be used to reject low hadronic energy events which 



are tagged through the vr — > /i — > e decay chain (see sec. 3.3.1 ). 

In order to achieve the rates and carry out the measurements given in Section |3] 
NuSOnG consists of a 3.5 kton (3 kton fiducial volume) isoscalar target with high seg- 
mentation resulting in good separation between electromagnetic and hadronic showers 
and muon tracks with good energy resolution for each. Good separation between 
hadronic and electromagnetic showers and good muon identification are necessary for 
separation of neutral and charged current events, and for low systematic errors on the 
measurements of the neutrino and antineutrino electron scattering cross sections. Fi- 
nally, good muon identification is critical for detecting inverse muon decay events for 
a precise flux measurement. 

Given the large size of the detector, ease of construction and low cost technologies 
are important. The long running time requires high stability and robust operation. 
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EM 
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Hadronic 
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Absorber 
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1.04//E 


0.72//E 


8% 
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0.80//E 


- 


5% 


2.8/8.3 Xo 
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CHARM II 


0.52/v^+0.02 


0.24//E+ 0.34 


5% 


0.5 Xo 
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digital) 












NuTeV 


0.86/\/^+0.022 


0.5/v^ +0.042 


10% 


5.8 Xo 


steel 


(Scintillator) 













Table 4: Comparison of high energy neutrino detectors. 



Our first design is shown in Figs. 14 and 15 and summarized in Table |3j NuSOnG 



consists of four calorimeters each with a muon spectrometer. 15 m decay volumes 
separate the four detector elements. Interspersing the decay volumes between the 
detectors will allow a calibration beam to be brought to each of the four detector 
regions. 

Each calorimeter has 500 Si02 2.5 cm (Xo/4) glass target planes interleaved with 
active detectors with two dimensional readout. The active detectors could be propor- 
tional tubes, scintillator panels, or a combination of both. These three options are 
discussed below. Neutrinos interact in the target planes, creating secondary particles; 
the active detector determines the total energies of the hadronic and electromagnetic 
secondaries. The muon detector measures the momentum of muon secondaries and 
serves to identify them. The pattern of the shower serves to identify the shower type: 
showers in which all the energy resides in ten of fifteen planes will be electromagnetic, 
and more extended showers will be hadronic. The lateral extent of the shower also 
resolves electromagnetic from hadronic showers. 

We have chosen an Si02 target. This material provides a balance between longer 
radiation length, important to particle ID issues, and shorter detector length, important 
for acceptance and calibration issues. The target could be commercial glass or thin 
walled plastic boxes filled with sand. Glass planes have the advantage of being easy 
to install and require no construction. Sand-filled boxes could be much less expensive. 
We will investigate both possibilities. Either way, Si02 has the advantage of being 
isoscalar ((iV„) > /(TV^) =0.998). Si02 has a density of 2.2 g/cm^; a high energy muon 
will lose 10 MeV per plane, which gives 5 GeV across all 500 planes in one calorimeter. 
Energy loss will also occur through electromagnetic showers. An example straight- 



through muon event from our initial GEANT4 detector simulation is shown in Fig. 16 
A michel electron with 30 MeV energy should be clearly visible across three planes. 
Each calorimeter is followed by a toroidal muon spectrometer consisting of magnetized 
iron plates interleaved with drift chambers. 

Other target materials, including emulsion, are under consideration, as has been 
discussed in previous sections of this document. These materials are not yet incor- 



49 



50 plane module 



24910.0 



4006.7 



Figure 14: NuSOnG detector showing calorimeter modules and muon detector. 
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Figure 15: The full NuSOnG experiment showing four detectors separated by a decay volume. 



Figure 16: A 100 GeV muon traversing the detector from the NuSOnG initial GEANT4 
Monte Carlo. 
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porated into the preliminary design presented here, but should be straightforward to 
include in the future. 

The design must address beam correlated backgrounds. These include backgrounds 
arising from debris (muons, remnants of hadronic showers) from neutrino interactions 
in the earth surrounding the detector. We plan for a forward veto consisting of a 
three layer scintillator hodoscope. Since our detector is so long, we may also need a 
veto system along the sides, top, and bottom of the calorimeters. We plan a Monte 
Carlo study of veto requirements in the coming months. Cosmic rays muons and 
their attendant showers present a beam-uncorrelated background which we will need 
to eliminate. We envisage a counter on the top of the detector similar to that used by 
the MINOS experiment. 

NuTeV showed the value of continuous beam calibration and this will be discussed 
in a separate section. 

While our detector is quite large, the robust, simple design will make the cost 
and construction manageable. The modules design makes upgrades and improvements 
straightforward. While we are designing with an initial four to five year run in mind, 
this detector can be put to other uses should the physics warrant. 



4.3.1 Active detector options 

The active detector performs two roles: first, it tracks the particles emerging from a 
neutrino interaction; second, it samples the particle's energy loss along the trajectory 
giving an measurement of the total energy. Simplicity, robustness, and high efficiency 
are essential, as is low cost. 

Two technologies immediately present themselves: gas-filled proportional tubes and 
plastic scintillator read out by phototubes. Both have been used in several experiments 
(see Table |4]). At this point, it is not clear to us which is the best approach for NuSOnG. 
We are also considering a design with both proportional tubes and scintillator. In the 
coming months, we plan to study the performance of each via simulation, develop 
preliminary design prototypes, and carry out a detailed cost estimate. We describe 
each detector concept below. 



4.3.1.1 Proportional Tubes 



A first design for a proportional tube active detector is shown in Fig. 17 Each active 
detector plane is made from five 1 m x 5 m extruded aluminum panels. Each panel 
contains fifty 1 cm x 2 cm drift cells. A 50 fim wire is strung down the center of each 
tube, and the applied high voltage produces both drift and proportional amplification 
fields. Ar:C02 (80:20) provides a good candidate for a fill gas; with 1.8-2 kV applied 
to the wire, the drift field will give a drift velocity of about 50 /im/ns and a gain of 
3000. A minimum ionizing particle crossing the 1 cm cell will deposit 2.7 keV of energy, 
liberating about 160 drift electrons in ten or so clusters. The drift time across the cell 
will be about 500 ns and proportional multiplication will give a collected charge of 80 
fC over a time of 250 ns. 

As an example of a readout scheme, we look to the ATLAS Transition Radiation 
Tracker (TRT) ASIC chips. The TRT readout has a peaking time of 7.5 ns and a charge 
threshold of 2 fC, making them well matched to our proportional tubes. Each chip 
set reads out sixteen channels and can be configured to provide trigger information. 
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Figure 17: Proportional tubes design. 



The TRT system, which is based around 6 mm straw tubes, has achieved a spatial 
resolution of 127 ^m, albeit with higher energy deposition resulting from the use of a 
xenon mixture. Scaling by the energy deposition gives a resolution of 200/im for our 
argon-filled tubes. The TRT readout chip set has sufficient charge sensitivity to allow 
us to use charge division; this should give position resolution of 5-10 cm along the wire. 



4.3.1.2 Scintillating Strips 

The second option uses planes of scintillator strips read out with green wavelength- 
shifting fibers fed into multi- anode photomultipliers. This option would be similar to 
that used for the SciBar detector in K2K, the Minos neutrino detector, and the Opera 
neutrino detector. NuSOnG would have 2500 5 m by 5 m planes with each plane made 
up of 128 3.9 cm x 1.3 cm strips. Each 64 strip plane will be separately wrapped in 
an Al skin that will provide the light seal and strength for the module. 

The scintillator strips will be coextruded with a Ti02 reflective coating and have a 
1.8 mm diameter hole in the middle. A 1.5 mm diameter green wavelength-shifting fiber 
will be put in the hole and routed to multianode photomultipliers for readout. The 64 
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wavelength-shifting fibers on one side of a plane will be coupled to a Hamamatsu M64 
multianode photomultiplier tube. The readout side will alternate between subsequent 
planes to improve uniformity. The fiber end opposite to the tubes will be polished and 
mirrored to increase the light output and uniformity. Planes will alternate between 
horizontal and vertical strips to provide two view tracking; readout tubes will alternate. 

The readout would be based on a custom ASIC combined with a standard FPGA. 
One example is the 64 channel MAR0C2 custom integrated circuit, designed at LAL 
(Orsay) for the ATLAS luminosity monitor. This chip allows adjustment of the elec- 
tronic gain of each of the 64 channels, which will be needed to correct for the expected 
factor of 3 pixel-to-pixel gain variation of the M64 tubes. The system provides a self- 
triggering analog readout into an external flash ADC. A fast discriminator signal for 
triggering is also available for each strip with a common threshold. 

Based on the performance of the SciBar detector, a minimum ionizing particle 
traversing a strip will yields about 20 photoelectons close to the tube, and the strip/fiber 
system will have an attenuation length of 3.5 m. This would then produce about 10 
photoelectrons at the center of the detector per plane. 

4.3.1.3 Hybrid Design 

Our initial estimates indicate the scintillator option may cost more than the pro- 
portional tube option. However, the scintillator system described above does provide a 
stable, easy to characterize active detector. In particular, scintillating strips offer very 
stable response that does not vary with pressure or temperature. We will investigate 
a hybrid system in which every fourth or eighth plane (one or two radiation lengths) 
would be a scintillator panel. The high granularity of the proportional tube design 
would give good pattern recognition, and the excellent energy resolution of the scin- 
tillator would give a better energy measurement. Reducing the fiber spacing in the 
scintillator may be possible; this would reduce the cost. 

One issue with adding 12-25% scintillator would be the change in the fraction of 
protons in the detector. The precise change depends on the scintillator used, but for 
CH4 and one scintillator panel every quarter radiation length, the proton-neutron ratio 
changes from 0.998 to 0.940. The impact of this change will have to be balanced against 
the cost reduction and stability improvement. This will be part of our Monte Carlo 
effort in the coming months. 

4.3.2 Toroid Spectrometers 

High energy muons produced in charged-current interactions will be momentum an- 
alyzed in three iron toroid spectrometers downstream of each subdetector (set of ten 
"stacks"). Each spectrometer will be composed of layers of magnetized iron instru- 
mented with drift chambers for tracking. 

Since NuSOnG will see muons of the same energies as NuTeV/CCFR a similar 
arrangement for measuring muon momenta would be suitable. CCFR used sections 
of 8" thick steel washers instrumented with scintillator hodoscopes for calorimeter]^ 

^The CCFR arrangement used two C-shaped sections with a horizontal crack at the center to allow 
placement of hall probes for field calibration. This crack would be eliminated in NuSOnG and instead small 
slots could be included for this purpose. 
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drift chamber stations drift chamber stations 



Figure 18: Conceptual Schematic for a NuSOnG toroid element. The three sections contain 
the same amount of steel (eight washers of 8" each). The upstream most section has addi- 
tional drift chamber stations to improve acceptance for low energy muons. Each of the five 
drift chamber stations has 3x and 3y view chambers. 

Tracking was performed using four views of each x and y chambers (0.5 mm coordinate 
resolution) in three gaps located after each 1.6 m of steel. The magnetic field was 
produced by four coils carrying approximately 1500 A each which passed through the 
center hole. The field was nearly radially symmetric and pointed in the azimuthal 
direction with magnitude ranging from 1.9 T near the center hole to 1.55 T near the 
outer edge (at R=1.8 m). Details can be found in reference |142j . 

Figure [TS] shows a possible arrangement for a NuSOnG toroid spectrometer. One 
"Upstream section" and two downstream "Standard sections" are shown. The down- 
stream sections contain eight 8" washers with one drift chamber station with 3x and 
3y view chambers each. The most upstream section of a spectrometer unit has two 
additional drift chamber stations to improve acceptance for low energy muons. To pass 
the coil through this arrangement the upstream chamber stations would be half size 
(the same chambers but rotated for each view). Each of the three sections contain the 
same amount of steel. Hodoscope paddles could be added in each chamber station for 
triggering purposes. Resolution of this arrangement would be dominated by multiple 
Coulomb scattering and would be ~ 11% independent of momentum. 

The NuSOnG arrangement will provide good acceptance for high energy primary 
muons of both signs since in a sign-selected beam the can be routinely operated with the 
polarity set to focus the primary muon. Very high energy particles can be tracked into 
the downstream target sections with a long lever arm and their momentum analyzed, 
(resolution for very high energy muons (> 150GeV was limited in NuTeV and CCFR; 
this resulted in large uncertainties in measuring flux in the high energy tail of the 
beam). Improving flux measurements in this region may help constrain kaon fluxes 
and therefore electron neutrino beam contamination. 
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4.3.3 Detector Calibration 



A thorough and precise cahbration of the entire detector will be required to achieve 
the physics goals of NuSOnG. Some of the response features of the detector can be 
understood using beam and cosmic ray muon samples, but a dedicated calibration effort 
will be required to study the hadronic and electromagnetic response of the detector 
and to measure the absolute energy scales. Precise calibration of a detector of this size 
will require a dedicated in situ calibration beam such as was used in NuTeV for this 
purpose [1431. 

The requirements for NuSOnG calibration beam would be similar to those of 
NuTeV. Tagged beams of hadrons, electrons, and muons over a wide energy range 
(5-200 GeV) would be required. The calibration beam should have the ability to be 
steered over the transverse face of the detector in order to map the magnetic field of 
each toroid with muons. This could be accomplished in several ways; for example, gaps 
of a few meters in front of each toroid could be incorporated into the design, and the 
beam could be steered into each toroid in turn; or the toroids could each be moved 
into the test beam for these calibration runs. Steering for hadrons and electrons would 
be less crucial than it was in NuTeV's case but would still be useful. 

The calibration beam can be constructed with a similar design to NuTeV. Upstream 
elements were used to select hadrons, electrons, or muons. An enhanced beam of 
electrons was produced by introducing a thin lead radiator into the beam and detuning 
the portion of the beam downstream of the radiator. A radiator was also used in the 
nominal beam tune to remove electrons. Particle ID (a threshold cerenkov and TRDs) 
was incorporated in the spectrometer and used to tag electrons when running at low 
energy. A pure muon beam was produced by introducing a 7 m long beryllium filter 
in the beam as an absorber. 

The NuTeV calibration spectrometer was able to determine incoming particle mo- 
menta with a precision of better than 0.3% absolute. This was accomplished by two 
means. First, precisely calibrated dipole spectrometer magnets were used, with J Bdi 
known to better than 0.1% in the region traversed by the beam. Secondly, the bend 
angle was determined to better than 0.1% using drift chambers positioned over the 
150 m spectrometer. This long lever arm allowed a modest alignment uncertainty of 
a few mm to translate into only a 0.1% uncertainty in the absolute momentum scale. 
The event-by-event resolution of the spectrometer, dominated by multiple scattering in 
the drift chamber walls, was better than 0.3% for most energies. (Helium in the region 
between the last dipole and the upstream part of the detector reduced the scattering 
in air). 

Figure 19 shows the NuTeV calibration beam configuration and the long lever arm 
spectrometer used to tag particle momenta with an absolute precision of better than 
0.3%. The most downstream dipole was mounted on a rotating stand which gave the 
ability to steer the beam out of the plane. 

The NuSOnG goal of the calibration precision would be to measure energy scales 
to a precision of about 0.5%. NuTeV achieved 0.43% precision on absolute hadronic 
energy scale and 0.7% on absolute muon energy scale (dominated by the ability to 
accurately determine the toroid map) . Precise knowledge of the muon energy scale is 
especially important in order to achieve high measurement accuracy on the neutrino 
fluxes using the low-u method. For example a 0.5% precision on muon energy scale 
translates into about a 1% precision on the flux. Both energy scales are important 
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Figure 19: (Top) Components of the NTEST beamline used to calibrate the NuTeV detector. 
Four different thicknesses of converter material at NTACON were used to select pure hadrons 
or electrons. The 7 m long Be filter(NTBBE) was used to select pure muons. The numbers 
on the left-hand-side of each component indicate the relative distance of the component to 
the primary target (NT8TGT) in meters. (Bottom) NuTeV's long lever arm spectrometer. 
The four dipole bend magnets were located in an enclosure approximately 70 m upstream of 
the Lab E detector. The spectrometer spanned over 150 m in length; this allowed precision 
measurement of the bending angle. 
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for precision structure function measurements and were the largest contributions to 
structure function measurement uncertainties in NuTeV [89j. 

4.4 Possible Locations 

Fig.|20|shows a possible location for the NuSOnG beam target hall and detector. Other 
layouts are possible; this is just meant to provide an example. 

This layout assumes that the beam is extracted at AO from the TeVatron and 
directed through the Switchyard Complex to a new targeting hall. This location allows 
low luminosity extraction down existing beamlines for the calibration beam. 

The detector is located near the New Muon Lab. This is a region with more than 
200 m of clear length, with roads and utilities nearby. 

The calibration beam could be delivered to the NuSOnG hall using a scheme similar 
to that used in NuTeV with the NTest beamline. Only a short extension of the existing 
NTest line would be required to reach a detector located near the New Muon Lab. The 
beam was split off from the same beamline (Ncenter) and then bent around to impinge 
on the detector at a 43 mrad angle. 
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Figure 20: Aerial view of Fermilab showing the Tevatron, external beam lines and potential 
site for NuSOnG target and detector halls. 
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5 Summary 



NuSOnG is an experimental program with high discovery potential. The precision neu- 
trino scattering measurements probe terascale physics and will complement discoveries 
at the LHC. Through precision electroweak measurements, NuSOnG will be sensitive 
to such new phenomena as extra Z bosons with masses beyond the 1 TeV (depending 
on the model) and compositeness scales above 5 TeV. The NuSOnG measurement of 
the coupling to the when combined with the LEP measurement of the invisible 
width, is a more sensitive method to search for new physics than this same measure- 
ment at the ILC. NuSOnG can also probe the existence of ncutrissimos, moderately 
heavy neutral heavy leptons which may be produced at the LHC, but which could be 
difficult to reconstruct and identify. A wide range of direct searches for new particles 
and interactions can be accomplished. The high neutrino flux and isoscalar target will 
make allow measurements which probe deeper into nuclear structure. 

The high energy neutrino facility, which uses 800 GeV protons from the TeVatron, 
has been endorsed by the 2007 Fermilab Steering Group. While NuSOnG is the first 
to propose an experiment for this facility, a wide range of interesting measurements 
can be made on this line. 

The proposed 3 kton (fiducial) NuSOnG detector design, which is opitmized for the 
physics goals, is based largely on the experiences of NuTeV and CHARM H. The basic 
technology is straightforward, although challenges exist because of the high precision 
demanded by the physics goals. Detailed simulations of the detector are now underway 

Our plan is to develop these ideas over the coming months. We plan to submit a 
proposal to the Fermilab Directorate in the near future. 
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